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Abstract 

We compute the leading-color (planar) three-loop four-point amplitude of N = 4 supersymmet- 
ric Yang-Mills theory in 4 — 2e dimensions, as a Laurent expansion about e = including the 
finite terms. The amplitude was constructed previously via the unitarity method, in terms of two 
Feynman loop integrals, one of which has been evaluated already. Here we use the Mellin-Barnes 
integration technique to evaluate the Laurent expansion of the second integral. Strikingly, the 
amplitude is expressible, through the finite terms, in terms of the corresponding one- and two-loop 
amplitudes, which provides strong evidence for a previous conjecture that higher-loop planar N = 4 
amplitudes have an iterative structure. The infrared singularities of the amplitude agree with the 
predictions of Sterman and Tejeda-Yeomans based on resummation. Based on the four-point result 
and the exponentiation of infrared singularities, we give an exponentiated ansatz for the maximally 
helicity- violating n-point amplitudes to all loop orders. The 1/e 2 pole in the four-point amplitude 
determines the soft, or cusp, anomalous dimension at three loops in N = 4 supersymmetric Yang- 
Mills theory. The result confirms a prediction by Kotikov, Lipatov, Onishchenko and Velizhanin, 
which utilizes the leading-twist anomalous dimensions in QCD computed by Moch, Vermaseren 
and Vogt. Following similar logic, we are able to predict a term in the three-loop quark and gluon 
form factors in QCD. 
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I. INTRODUCTION 



Maximally supersymmetric N = 4 Yang-Mills theory (MSYM) in four dimensions has a 
number of remarkable properties. There are good reasons to believe that, in the 't Hooft 
(planar) limit of a large number of colors N c , higher-loop orders are surprisingly simple [1]. 
In particular, the anti-de Sitter/conformal field theory (AdS/CFT) correspondence suggests 
a simplicity in the perturbative expansion of planar MSYM as the number of loops in- 
creases [1]. The Maldacena conjecture [2] states that the planar limit of MSYM at strong 
coupling is dual to weakly-coupled gravity in five-dimensional anti-de Sitter space. Based 
on this conjecture, one might expect observables in the strongly-coupled limit of MSYM to 
have a relatively simple form, due to the interpretation in terms of weakly-coupled gravity. 
On the other hand, the strong-coupling limit of a typical observable receives contributions 
from infinitely many terms in the perturbative expansion, as well as non-perturbative effects. 
How might the perturbative series be organized to produce a simple strong-coupling result? 
Some quantities are protected by supersymmetry — non-renormalization theorems lead to 
zeroes in the perturbative series, which certainly can bring about this simplicity [3, 4]. It 
has been less clear how the perturbative series for unprotected quantities might have the 
required simplicity [4, 5, 6]. 

One suggestion, confirmed through two loops for dimensionally-regulated on-shell scatter- 
ing amplitudes, is that an iterative structure exists [1], which may eventually allow the per- 
turbative series to be resummed into a simple result. In particular, the planar four-point two- 
loop amplitude of MSYM was shown to be expressible in terms of the corresponding one-loop 
amplitude. Roughly speaking (see eq. (4.1) for the precise formula), the two-loop amplitude 
is given by the square of the one-loop amplitude, plus a term proportional to the one-loop 
amplitude evaluated in a slightly different dimension, plus a constant. This result was found 
using the two-loop integrand [7, 8] obtained via the unitarity method [9, 10, 11, 12, 13], and 
the Laurent expansion in e = (4 — d)/2 of the associated two-loop planar box integral [14]. 

On-shell loop amplitudes in massless gauge theory have severe infrared (IR) singularities, 
arising from soft and collinear loop momenta. Regulated dimensionally, the singularities 
produce poles in the limit e — > 0, beginning at (D(e~ 2L ) for an L-loop amplitude. The two- 
loop iterative relation holds from (9(e~ 4 ) through O(e ), but it does not hold at 0(e 1 ). This 
observation is consistent with intuition that a simple structure need only exist near four 
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dimensions [1], where MSYM is a conformal theory, and where it should be dual to a gravity 
theory in anti-de Sitter space. 

Splitting amplitudes are functions governing the behavior of scattering amplitudes as 
two momenta become collinear. The two-loop splitting amplitude in MSYM has an iterative 
structure very similar to that of the four-point amplitude [1, 13]. Based on this structure, 
an iterative ansatz for the planar n-point two-loop amplitudes can also be constructed. The 
ansatz is very likely to be true for the maximally helicity violating (MHV) amplitudes (those 
with two negative helicities and the rest positive) because it ensures that these amplitudes 
have the correct factorization behavior in all channels. (For non-MHV amplitudes one would 
also need to ensure that the structure of the multi-particle poles is correct.) 

Amplitudes for scattering of on-shell massless quanta have considerable practical rele- 
vance, in the applications of perturbative QCD to collider physics. At the perturbative 
level, MSYM is a close cousin of QCD, although its amplitudes have a much simpler an- 
alytic structure, allowing their computation typically to precede the QCD result. In fact, 
the surprisingly simple structure of MSYM loop amplitudes has been unfolding for quite a 
while, beginning with the superstring-based evaluation of the one-loop four-point amplitude 
by Green, Schwarz and Brink [15]. Compact results for the n-point MHV amplitudes [9], 
and for all helicity configurations at six points [10], were among the early applications of the 
unitarity method of Dunbar, Kosower, and two of the authors [9, 10, 11, 12, 13]. Because the 
unitarity method builds amplitudes at any loop order from on-shell lower-loop amplitudes, 
any simplicity uncovered at the tree and one-loop levels should induce a corresponding ad- 
ditional simplicity at higher-loop orders. Indeed, the simplicity observed in the multi-loop 
four-point MSYM loop integrands (prior to performing loop integrations) was found in this 
way [7, 8]. 

Witten has proposed a duality between MSYM and twistor string theory [16], generalizing 
Nair's earlier description [17] of MHV tree amplitudes. This proposal, and the investigations 
it has stimulated into the structure of tree [18] and one-loop [19, 20] gauge theory scattering 
amplitudes, provide additional strong support for the notion that amplitudes — particularly 
MSYM amplitudes — should be remarkably simple. 

These results, particularly the two-loop iterative relation, lead to the natural conjecture 
that an iterative structure should continue to hold for higher-loop planar MSYM ampli- 
tudes [1]. The purpose of this paper is to verify the conjecture at the level of the three-loop 
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four-point amplitude, and to flesh out more of the likely structure beyond three loops. 

The planar three-loop four-point MSYM amplitude was found in ref. [7] via the unitarity 
method, and expressed in terms of just two independent loop integrals. To check for an 
iterative relation, we must first compute the expansion of these two integrals around e = 0, 
from the most singular terms, C(e -6 ), through the finite terms, O(e ). Fortunately, there 
has been much progress in multi-loop integration over the past few years [14, 21, 22]. One of 
the two integrals we need, a three-loop ladder integral, was computed through finite order 
recently by one of the authors [23], using a multiple Mellin-Barnes (MB) representation. 
In this paper, we present the expansion of the single remaining integral — and thus the 
expansion of the three-loop amplitude — through the finite terms. 

We wish to compare this expression with the expansions of products of one- and two- 
loop amplitudes. For this purpose, we must expand the one- and two-loop amplitudes to 
(9(e 4 ) and 0(e 2 ) respectively, which is two higher orders in e than was necessary at two 
loops. All of the expansions are given in terms of harmonic polylogarithms [24, 25]. We 
use identities to reduce the harmonic polylogarithms to an independent basis set. Taking 
into account intricate cancellations between the different amplitude terms, we find that the 
planar three-loop four-point amplitude does indeed have a simple iterative structure (see 
eq. (4.4)). 

To guide us toward the correct iterative relation, we employed properties of the three-loop 
amplitude's IR singularities [26], which must be respected by any such relation. In general, 
the IR singularities of loop amplitudes in gauge theory can be represented in terms of 
universal operators, acting on the same scattering amplitudes evaluated at lower loop order, 
as was first discussed at one and two loops [27, 28]. These operators are related to the soft (or 
cusp) anomalous dimension and other quantities entering the Sudakov form factor [29, 30], as 
was clarified recently [26] . The latter quantities play an important role in the resummation 
and exponentiation of large logarithms near kinematic boundaries, such as the threshold 
(x — > 1) logarithms in deep inelastic scattering or the Drell-Yan process [30, 31, 32]. 

In other words, the IR divergence structure of loop amplitudes are a priori predictable, 
up to sets of numbers (e.g. soft anomalous dimensions) that must be obtained by specific 
computations. Our four-point computation simultaneously provides a verification of the 
three-loop IR divergence formula [26], and a direct determination of two of the numbers 
entering it, for planar MSYM: the three-loop coefficients of the soft anomalous dimension 
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and of the (/-function for the Sudakov form factor [26, 30]. 

The three-loop four-point iterative relation, combined with information about how IR 
singularities exponentiate [30], and the factorization properties used at two loops [1], leads 
us to an exponentiated ansatz for the planar n-point MHV amplitudes at L loops. This 
ansatz naturally produces each loop amplitude as an iteration of lower-loop amplitudes, up 
to a set of constants which are as yet undetermined beyond three loops. (Two rational 
numbers at three loops are also undetermined.) By taking collinear limits of the ansatz, we 
obtain, as a by-product, an iterative ansatz for the L-loop splitting amplitudes of MSYM. 

We use the universal form of the divergences to define IR-subtracted finite remainder am- 
plitudes. (Similar subtractions are made in perturbative QCD when constructing finite cross 
sections for infrared-safe observables.) For our exponentiated ansatz, the finite remainder at 
L loops is strikingly simple: it is a polynomial of degree L in the one-loop finite remainder. 
This result applies directly to the finite remainder of the three-loop four-point amplitude, 
for which it follows from actual computation, not an ansatz. 

Infrared singularities provide a link between the scattering amplitudes computed here 
and the anomalous dimensions of gauge-invariant composite operators in MSYM, studied 
in the context of the AdS/CFT correspondence [4, 5, 33, 34]. Specifically, at three loops, 
the coefficient of the 1/e 2 IR singularity is controlled by the high-spin, or soft, limit of the 
leading- twist anomalous dimensions [26]. Equivalently, it appears in the x — > 1 limit of 
the kernels for evolving parton distributions fi(x,Q 2 ) in the scale Q 2 . The x — > 1 limit 
of the splitting kernels corresponds to multiple soft gluon emission, and is related to the 
soft (or cusp) anomalous dimension associated with a Wilson line [35]. The three-loop soft 
anomalous dimension in QCD has been computed by Moch, Vermaseren and Vogt as part 
of the heroic computation of the full leading-twist anomalous dimensions [36]. (The terms 
proportional to Nf were computed earlier [37].) 

The QCD result has been carried over to MSYM by Kotikov, Lipatov, Onishchenko 
and Velizhanin (KLOV) [38], using an inspired observation that the MSYM results may be 
obtained from the "leading-transcendentality" contributions of QCD. For the soft anomalous 
dimensions, which are polynomials in the Riemann ( values, ( n = ((n), the degree of 
transcendent ality is tallied by assigning the degree n to each ( n . The KLOV observation 
applies to the anomalous dimensions for any spin j; a similar accounting of harmonic sums 
Sm(j) is used to assign transcendentality in that case. Very interestingly, the three-loop 
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MSYM anomalous dimensions of KLOV were confirmed by Staudacher [39] through spin 
j = 8, building on earlier work of Beisert, Kristjansen and Staudacher [34] at j = 4, by 
assuming integrability and using a Bethe ansatz. Our determination of the three-loop soft 
anomalous dimension in MSYM now provides an independent confirmation of the KLOV 
result in the limit j — > oo. 

The iterative structure of MSYM is presumably tied to the issue of integrability of the 
theory [33, 34]. There has also been an interesting hint of a similar structure developing 
in the correlation functions of gauge-invariant composite operators in MSYM [40]; but its 
precise structure, if it exists in this case, has not yet been clarified. 

This paper is organized as follows. In section II we review known results for planar 
loop amplitudes in MSYM, focusing on the construction of the three-loop integrand for 
the four-point amplitude. The methods used to evaluate the two three-loop integrals are 
described in section III. In section IV we describe the iterative relation for the three-loop 
four-point amplitude. Then we present an exponentiated ansatz which extends the relation 
to n-point MHV amplitudes at an arbitrary number of loops. We discuss the consistency 
of this ansatz with exponentiation of infrared singularities. The consistency of our ansatz 
under factorization onto kinematic poles, particularly the collinear limits, is discussed in 
section V. In section VI we relate the anomalous dimensions and Sudakov coefficients 
appearing in the L-loop amplitudes to previous work in QCD and MSYM. Our conclusions 
are given in section VII. Appendix A summarizes properties of harmonic polylogarithms, 
while appendix B contains the results for all loop integrals encountered in our calculation 
of the amplitudes. 

II. GENERAL STRUCTURE OF MSYM LOOP AMPLITUDES 

It is convenient to first color decompose the amplitudes [12, 41] in order to separate 
the color from the kinematics. In this paper we will discuss only the leading-color planar 
contributions. These terms have the same color decomposition as tree amplitudes, up to 
overall factors of the number of colors, N c . The leading-A c contributions to the L-loop 
SU (N c ) gauge-theory n-point amplitudes may be written in the color-decomposed form as, 

L 
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where 7 is Euler's constant, and the sum runs over non-cyclic permutations of the external 
legs. In this expression we have suppressed the (all-outgoing) momenta hi and helicities Aj, 
leaving only the index i as a label. This decomposition holds for all particles in the gauge 
super- multiplet which are all in the adjoint representation. The advantage of this form is 
that the color-ordered partial amplitudes A n are independent of the color factors, cleanly 
separating color and kinematics. We will not discuss the subleading-color contributions here 
because there does not appear to be a simple iterative structure present for them [1]. 

In general, loop amplitudes in massless gauge theory, including MSYM, contain IR sin- 
gularities. This implies that a textbook definition of the S-matrix with fixed numbers of 
elementary particles does not exist. To define an S'-matrix in massless gauge theory, dimen- 
sional regularization — which explicitly breaks the conformal invariance — is commonly 
used. Once the universal IR singularities are subtracted, the four-dimensional limit of the 
remaining terms in the amplitudes may then be taken. In QCD, after combining real emis- 
sion and virtual contributions, these finite remainders are the quantities entering into the 
computation of infrared-safe physical observables [42]. It is worth noting that the finite 
remainders should also be related to perturbative scattering matrix elements for appropri- 
ate coherent states (see e.g. ref. [43]). The IR singularities for MSYM that we discuss in 
this paper are closely connected to those of QCD and are, in fact, a subset of the QCD 
divergences. As is typical in perturbative QCD, the S'-matrix under discussion here is not 
the one for the true asymptotic states of the four-dimensional theory, but for elementary 
partons. 

The unitarity method [9, 10, 11, 12, 13] provides an efficient means to obtain the inte- 
grands needed for constructing loop amplitudes. In this approach, the integrands for loop 
amplitudes are obtained directly from on-shell tree amplitudes without resorting to an off- 
shell formalism. A key advantage is that the building blocks used to obtain the amplitudes 
are gauge invariant and posses simple properties under extended supersymmetry, unlike 
Feynman diagrams. (Implicit in this approach is the use of a supersymmetric regulator, 
such as the four-dimensional helicity (FDH) scheme [44], a variation on dimensional reduc- 
tion (DR) [45].) The unitarity method derives its efficiency from the ability to use simplified 
forms of tree amplitudes to produce simplified loop integrands. 

The unitarity method expresses the amplitude in terms of a set of loop integrals. Experi- 
ence shows that such integrals can be evaluated in terms of generalized polylogarithms. At 
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FIG. 1: The result for the leading-color two-loop amplitude in terms of scalar integral functions, 
given in eq. (2.4). 

one loop a complete basis of dimensionally regularized integral functions is known [9, 10, 46], 
in general, reducing the integration problem to that of determining coefficients of the basis 
integrals. For four-point amplitudes only a single scalar box integral appears. At two and 
higher loops an analogous basis of integral functions is not known, and the integrals must 
be evaluated case by case. The two-loop massless planar double-box integral has, however, 
been evaluated in ref. [14] and is given in terms of harmonic polylogarithms [24, 25] through 
0(e 2 ) in eq. (B6) of the second appendix. One of the integrals appearing in the three-loop 
four-point amplitude has also been previously evaluated [23], and is given in eq. (B8). 

The one-loop four-point amplitude in MSYM was first calculated by taking the low energy 
limit of a superstring [15]. After scaling out a factor of the tree amplitude via, 



Here l\ is the one-loop scalar box integral, multiplied by a convenient normalization fac- 
tor, and is defined in eq. (Bl) of appendix B I. This box integral is identical to the one 
encountered in scalar 3 theory. Its explicit value in terms of harmonic polylogarithms is 
given through C(e 4 ) in eq. (B3). We keep the higher-order terms in e because they will 
contribute when we write the three-loop amplitude in terms of the one- and two-loop ampli- 
tudes. The factor of 1/2 in eq. (2.3) is due to our normalization convention for exposed 
in eq. (2.1) where a compensatory "2" appears in the brackets. This convention follows the 
QCD literature on two-loop scattering amplitudes (see e.g., ref. [47]). 



(2.2) 



the result for the one-loop four-point amplitude is rather simple, 



M i l\e) = - l -stl ( i\s,t). 
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FIG. 2: The rung insertion rule for generating higher- loop integrands from lower loop ones, given 
in ref. [7]. 

The two-loop MSYM four-point amplitudes were obtained in ref. [7] using the unitarity 
method, with the result for the planar contribution, 



which is schematically depicted in fig. 1. The two- loop scalar integral I\ is defined in 
eq. (B4). The scalar double box integral I^\s, t) was first evaluated through O(e ) in terms 
of polylogarithms by one of the authors using multiple MB representations [14]. In eq. (B6), 
we give this integral through 0(e 2 ). The higher-order terms in e are again needed because 
they will appear in our iterative relation for the three-loop amplitude. The result (2.4) has 

been confirmed using the two-loop four-gluon QCD amplitude for helicities ( h+) [47], 

which can be converted into the four-gluon amplitude in MSYM by adjusting the number 
and color of states circulating in the loop [1]. 

The original calculation [7] of the coefficients of the integrals in eq. (2.4) used iterated 
two-particle cuts, which are known to be exact to all orders in e since they involve precisely 
the same algebra used to obtain the one-loop amplitude (2.3). Beyond two loops, an ansatz 
for the planar contributions to the integrands was proposed in terms of a "rung insertion 
rule" [7, 8]. This ansatz was based on an analysis of two- and three-particle cuts, as well as 
cuts with an arbitrary number of intermediate states, but where the intermediate helicities 
are restricted so that the amplitudes on either side of the cut are MHV amplitudes. At three 
loops, the planar integrals generated by the rung rule can be constructed using iterated two- 
particle cuts, so the ansatz is reasonably secure. However, beyond three loops (and even at 
three loops for non-planar contributions) the rung rule generates diagram structures that 
cannot be obtained using iterated two-particle cuts. It is less certain that the rung rule gives 
the correct results for such contributions. There are also potential contributions coming from 
(— 2e)-dimensional parts of loop momenta, which have been dropped in the analysis of the 
three-particle and MHV cuts. These contributions would need to be kept in order to prove 




(2.4) 
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rigorously that the rung rule correctly gives all contributions. 

It is worth noting that while the integrand obtained from the rung insertion rule is only 
an ansatz, the results of this paper provide strong evidence that it is the complete answer, 
at least for the planar contributions at three loops. As we shall discuss in section IV B, the 
IR divergences of eq. (2.5) are fully consistent with the known form of the three-loop IR 
divergences [26]. Moreover, the non-trivial cancellations required by the iterative relations 
described in section IV imply that there are no missing pieces. 

In any case, we use the rung rule as our starting point for evaluating the planar three-loop 
MSYM amplitudes. According to this rule one takes each diagram in the L-loop amplitude 
and generates all the possible (L + l)-loop diagrams by inserting a new leg between each 
possible pair of internal legs as shown in fig. 2. From this set the diagrams which have 
triangle or bubble subdiagrams are removed. The new loop momentum is integrated over, 
after including an additional factor of + / 2 ) 2 in the numerator, where l\ and l 2 are the 
momenta flowing through each of the legs to which the new line is joined, as indicated in 
fig. 2. Each distinct (L + l)-loop contribution should be counted once, even though it can 
be generated in multiple ways. (The contributions which correspond to identical graphs but 
have different numerator factors should be counted as distinct.) The (L + l)-loop planar 
amplitude is then the sum of all distinct (L + l)-loop diagrams. The diagrams generated by 
the iterated two-particle cuts have an amusing resemblance to Mondrian's artwork; hence it 
is natural to call them "Mondrian diagrams." 

Applying this rule to the three-loop planar amplitude gives the explicit form of the 
integrand [7], 

Mf(e) = ~st (s 2 li 3)a (s,t) + 2sli 3)b (t,s) +t 2 /f a (t, s ) + 2t if )6 ( S , t)) . (2.5) 

This integrand is depicted in fig. 3 [70]. The second and third integrals in the figure are 
equal, as are the fifth and sixth, accounting for the appearance of six diagrams in fig. 3, but 
only four terms in eq. (2.5). The integrals lf^ a and I^ b appearing in the amplitude are 
defined in eqs. (3.1) and (3.2). The first of these integrals has been evaluated in ref. [23]. The 
evaluation of the second integral is outlined in the next section. The expansions of these 
integrals through O(e ), in terms of harmonic polylogarithms, are presented in eqs. (B8) 
and (BIO). 
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FIG. 3: Mondrian diagrams for the three-loop four-point MSYM planar amplitude given in eq. (2.5). 
The second and third diagrams have identical values, as do the fifth and sixth. The factors of 
(7i + I2) 2 denote numerator factors appearing in the integrals, where l\ and I2 are the momenta 
carried by the lines marked by arrows. 
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FIG. 4: The two integrals appearing in the three-loop amplitude. The "ladder" integral (a) has 
no factors in the numerator. The "tennis court" integral (b) contains a factor of (p + r) 2 in the 
numerator. 

III. EVALUATING TRIPLE BOXES 



The two three-loop integrals appearing in the four-point amplitude (2.4), and depicted 
in fig. 4 are 

d d pd d rd d q 



I ( l )a (s,t) = (-te^n' d / 2 f J 



p 2 {p — ki) 2 (p — ki — k 2 ) 2 
1 



x 



(p + r) 2 r 2 (q — r) 2 (r — k% — /c 4 ) 2 q 2 (q — A; 4 ) 2 (q — k% — /c 4 ) 
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and 



ii 3)b ( s ,t) = h 




d d pd d r d d q (p + r) 2 



p2 g2 r 2 (j) _ (j ) J rr _ ^2 

1 



(3.2) 
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(p + r - fci - k 2 ) 2 (p + r + k 4 ) 2 (q - k 4 ) 2 (r + p + q) 2 (p + q) 



2 ' 



where dimensional regularization with d = 4 — 2e is implied. 

The ladder integral, lj^ a , was evaluated in ref. [23], in a Laurent expansion in e up to the 
finite part, by means of the strategy based on the MB representation which was suggested in 
ref. [14] and applied for the evaluation of the massless on-shell double boxes. This strategy 
is presented in detail in Chapter 4 of ref. [48]. Here its basic features are briefly summarized. 

The strategy starts with the derivation of an appropriate multiple MB representation. 
MB integrations are introduced in order to replace a sum of terms raised to some power by 
their products raised to certain powers, at the cost of having extra integrations: 



where —Re A < (3 < 0. The simplest possible way of introducing an MB integration is 
to write down a massive propagator as a superposition of massless ones. In complicated 
situations, one starts from Feynman or alpha parameters and applies (3.3) to functions 
depending on these parameters. Of course, it is natural to try to introduce a minimal 
number of MB integrations. Anyway, after introducing sufficiently many MB integrations, 
one can evaluate all internal integrals over Feynman/alpha parameters in terms of gamma 
functions and arrive at a multiple MB representation with an integrand expressed in terms 
of gamma functions in the numerator and denominator. 

It turns out to be very convenient to derive a multiple MB representation for loop- 
momentum integrals of a given class with general powers of the Feynman propagators. Such 
a general derivation provides a lot of crucial checks and can then be used for any integral of 
the given class. Moreover, it provides unambiguous prescriptions for choosing contours in 
MB integrals, where the poles with T(. . .— z) dependence are to the right of the integration 
contour and the poles with F(. . . + z) dependence are to the left of it. 

To evaluate a given Feynman integral represented in terms of a multiple MB integral in 
an expansion in e one needs first to understand how poles in e are generated. A simple 
example is given by the product T(e + z)T(—z) which generates the singularity at e — > 
because, in this limit, there is no place for a contour to go between the first left and right 




(3.3) 
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poles of these two gamma functions, at z = — e and z = 0, respectively. To make the singular 
behavior in e manifest one can integrate instead over a new contour where the pole at z = — e 
is to the right of the contour (for example, (5 — — 1 in eq. (3.3), where A = e is assumed 
to have a small positive real part), plus a residue at this pole. We refer to the integral 
over the new contour as "changing the nature" of the first pole of T(e + z). In complicated 
situations, singularities in e are not visible at once, after one of the MB integrations. To 
reveal them one uses the general rule according to which the product T(a + z)T(b — z), with 
a and b depending on other MB integration variables, generates, due to integration over z, 
a singularity of the type F(a + b). 

Thus, to reveal the singularities in e one analyzes various products of gamma functions 
in the numerator of a given integrand, implying various orders of integration over given 
MB variables. After such an analysis, one distinguishes some key gamma functions which 
are responsible for the generation of poles in e. Then one begins the procedure of shifting 
contours and taking residues, starting from one of these key gamma functions. After taking 
a residue, one arrives at an integral with one integration less; one then performs an analysis 
of the generation of singularities in e in the same spirit as for the initial integral. For the 
integral with the shifted contour, one takes care of a second key gamma function in a similar 
way. As a result of this procedure, one obtains a family of integrals for which a Laurent 
expansion of the integrand is possible. To evaluate these integrals expanded in e, up to some 
order, one can use the second and the first Barnes lemma and their corollaries. A collection 
of relevant formulae are given in Appendix D of ref. [48]. 

The technique of multiple MB representation has turned out to be very successful, at least 
in the evaluation of four-point Feynman integrals with two or more loops and severe soft 
and collinear singularities (see refs. [14, 23, 49, 50]), so that it is natural to apply it to the 
evaluation of the three-loop tennis-court integral (3.2), which is the only missing ingredient 
of our calculation. Let us outline the main steps, following the strategy characterized above. 

An appropriate MB representation can be derived straightforwardly, in a way similar to 
the treatment of the ladder triple box integral (3.1) in ref. [23]. Indeed, one can derive an 
auxiliary MB representation for the double box with two legs off shell, apply it to the double 
box subintegral in (3.2), and then insert it into the well-known MB representation for the 
on-shell box (see, e.g., Chapter 4 of ref. [48]). As a result, an eightfold MB representation 
can be derived for the general diagram of fig. 4b with the eleventh index corresponding to 
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the numerator [{p + r) 2 ] 011 . For our integral with the powers a\ = . . . = a w = 1 and 
an = —1, this gives 

r(3)6 



r(-2e)(-s) 1 + 3e t 2 



r(-3e - w)Y{\ + Zl + z 2 + z 3 )T(-l -e-z l - z 3 )Y{\ + z 1 + z 4 ) 
T(l - z 2 )T(l - z 3 )T(l - z 6 )T(l -2e + Zl + z 2 + z 3 ) 
r(— 1 - e - z\ — z 2 — z 4 )Y(2 + e + z\ + z 2 + z 3 + z 4 ) 
X r(-l - 4e - ^ 5 )r(l - 24 - ^ 7 )r(2 + 2e + z 4 + z 5 + z 6 + z 7 ) 
xT(-e + z l + z 3 - z 5 )Y(2 — w + z 5 )T(-l + w - z 5 - z 6 ) 

xT(z 5 + z 7 - 2i)r(i + 2; 5 + ^ 6 )r(-i + w - 24 - -25 - z 7 ) 

xT(-e + zi + z 2 - z 5 - z 6 - z 7 )T(l - e - w + z 4 + z 5 + z 6 + z 7 ) 

xT(l + e - zi - z 2 - z 3 + z 5 + z 6 + z 7 ) . (3.4) 

There is a factor of T(— 2e) in the denominator, so that the integral is effectively sevenfold. 

A preliminary analysis shows that the following two gamma functions are crucial for the 
generation of poles in e: 

r(— 1 + w - z 5 - 2 6 )r(— 1 + w - z 4 - z 5 - z 7 ) . (3.5) 

The first decomposition of (3.4) reduces to taking residues and shifting contours with respect 
to the first poles of these two functions. We obtain 

T = lf )h = Too + T i + T10 + Tn ■ (3.6) 

The term T i denotes minus the residue at z 7 = — 1 + w — z 4 — z^ and changing the nature of 
the first pole of T(— — z 5 — zq); the term Ti denotes minus the residue at z e = — 1+w— z 5 
and changing the nature of the first pole of r(— 1+w— z 4 — z^— z 7 ); the term Tn corresponds to 
taking both residues; and T o refers to changing the nature of both poles under consideration. 

For each of these four terms, one proceeds further using the strategy of shifting contours 
and taking residues. One can arrive at contributions which are labelled by sequences of 
gamma functions. Let us denote by f (. . . ± taking the residue at the first pole of this 
gamma function with respect to the variable z iy and by T*(. . . ± changing the nature 
of this pole. If T(. . . ± Zj) participates then both variants are implied. If there is only one 
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^-variable in an argument of a gamma function then it is not underlined. The contributions 
that start from order e 1 in the Laurent expansion are not listed. So, for Too, one can arrive 
at the following eleven contributions: 

{f (-1 -e-z x - z 3 ), r(-e + z 2 ), f (-1 - 2e - z 5 ), F(-2e + z 6 ), f (-e + z 3 + z 7 ), 
f(-e + z 3 ),r(-2e + z 4 )}, 

{f (-1 - e - z, - z 3 ),T(-e + z 2 ),T{-l - 2e - z 5 ), f (-2e + z 6 ), f (-e + z 3 + z 7 ), 
T*(-e + z 3 ),T(-e-z 3 + z A )}, 

{f (-1 - e - z x - z 3 ), f (-e + z 2 ), f (-1 - 2e - * 5 ), f (-2e + z 6 ), 
r*(-e + 2 3 + z 7 ), f (-e - 2 3 + ^)}, 

{f (-1 - e - z x - z 3 ), f (-e + z 2 ), f (-1 - 2e - z 5 ), r*(-2e + z 6 ), f (-e + ^ 3 + z 7 ), 
f(-e + z 3 ),r(-2e + z 4 )}, 

{f (-1 -e-z x - z 3 ), f (-e + z 2 ), - 2e - z 5 ), f (-e - z 3 + zj}, 

{f (-1 - e - a - z 3 ), T*(-e + z 2 ), f (-1 - 2e - z 5 ), f (-2e + z 6 ), f (-e + z 3 + z 7 ), 

f(-e + 2; 3 ),r(-2e + Z4)}, 

{f (-1 - e - a - * 3 ), r*(-e + z 2 ), f (-1 - 2e - z 5 ), T*(-2e + z 6 ), f (-e + z 3 + z 7 ), 
f(-e + z 3 ),f(-2e + z 4 )}, 

{r*(-i - e - ^ - z 3 ), r(-e + ^ + 23 - z 5 ), r(-e + 23 + z 7 ), r(-e + * 3 )}. 

The rest of the 203 contributions present in T i + T 10 + T u can be described in a similar 

way. 

The final result for (3.2) is presented in eq. (BIO) of appendix B3. The evaluation of this 
integral has turned out to be rather intricate. The level of complexity is roughly five times 
the corresponding complexity of the ladder triple box. Therefore, systematic checks are quite 
desirable. A powerful independent check can be provided by evaluating the leading orders 
of the asymptotic behavior in some limit. Indeed, such checks were essential in previous 
calculations — see refs. [14, 23, 50]. Here we shall outline an independent evaluation of the 
dominant terms in eq. (BIO) in the limit s/t — > 0. 

The limit s/t — > is of the Regge type which is typical of Minkowski space. Hence the 
well-known prescriptions for limits typical of Euclidean space, written in terms of a sum 
over subgraphs of a certain class (see refs. [51, 52]), are not applicable here. However, one 
can use more general prescriptions formulated in terms of the so-called strategy of expansion 
by regions [52, 53, 54]. This approach is universal and applicable for expanding any given 
Feynman integral in any asymptotic regime. 
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An essential point of this strategy is to reveal regions in the space of the loop momenta 
which generate non-zero contributions. A given region is characterized by some relations 
between components of the loop momenta. In particular, in the case of our limit s/t — > 0, 
in the region where all the loop momenta are hard, all the components of the loop momenta 
are of order \ft. It turns out that the most typical regions relevant to the Regge and 
Sudakov limits are 1-collinear (lc) and 2-collinear (2c) regions. (Here "a-collinear" means 
that an appropriate loop momentum is collinear with external leg a.) The crucial part of 
the strategy of expansion by regions [53] is to expand the integrand in a Taylor series in 
parameters which are small in a given region and then extend the integration to the whole 
space of the loop momenta, i.e., forget about the initial region. Another prescription of 
this strategy is to put to zero any integral without scale (even if it is not regularized by 
dimensional regularization). 

In the case of the ladder triple box (3.1), in the Regge limit t/s — > 0, only the (lc-lc-lc) 
and (2c-2c-2c) regions participate in the leading power-law behavior [55]. 

For the tennis-court integral (3.2), the evaluation procedure outlined above is formulated 
in such a way that the leading terms of the expansion at s/t — > can be clearly distinguished. 
All these terms arise after taking residues with respect to the variable w at w = or w = e. 
It turns out that only one contribution to the result (BIO) arises after taking a residue at 
w — e. It involves no integration (i.e., it is obtained from (3.4) by taking consecutively eight 
residues), so that it can be expressed in terms of gamma functions for general values of e: 

I^ b —(s, t) = ( _ g)1+ f ( 7 _ t)2 ^ e r (- £ ) 3r ( g ) 2r ( 1 + 2 ^) 2 • (3-7) 

It turns out that this term is nothing but the (lc-4c-us) contribution within the expansion 
by regions. It is generated by the region where the momentum of the line between the 
external vertices with momenta A; 2 and k% is considered ultrasoft (us), the loop momentum 
of the left box subgraph is considered 1-collinear and the loop momentum of the right box 
subgraph is considered 4-collinear — see fig. 4(b). (Details of the expansion in the Sudakov 
and Regge limits within the expansion by regions can be found in ref. [54] and Chapter 8 of 
ref. [52].) 

The rest of the contributions to the leading power-law behavior in the limit s/t — > 
correspond to taking residues at w — 0. They can be identified as the sum of the (lc-lc- 
lc) and (4c-4c-4c) contributions. The (4c-4c-4c) contribution can be represented by the 
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following fivefold MB integral: 



r(3)6,4c-4c-4c 
4 



(S,t) 



r(l + 3e + a; 2 )r(-3e-x 2 ) 



x 



1 




(_ s )i+3e+* 2 (_t)2+xx r (_ 2e _ Xl )r(i + xi)r(i + x 2 ) 
+lco 5 r(1 + ^ )r(1 + e + ^ + 23)r( _ 22)r( _ e _ 23 ) r (_ 23 ) 



(2tt«) 




r(— i - e - 21 - ^2)r(i - xi + zi + 2 2 + 2 3 )r(— l + xi — ^ — z 4 ) 

r(l - 2e - xi - x 2 + + 2 2 + z 3 )r(— 2 — 4e — x 2 — ^i — -24) 
T(l + z 4 )T(e + x 2 - z 2 - z 3 - z 5 )r(-2 + e + x 2 - zi - z 2 - z 3 - Zi - z 5 ) 
r(l + e + rr 2 - z 2 - 2 3 - z 5 )r(l + e + z 2 + z 3 + z 5 ) 
xT(l + z 5 )r(-l - e - rr 2 + z 2 - z A )T(-2e - x 2 + z 2 + z 3 + z 5 ) 

xT(2 - e- Xl + zi + z 2 + z 3 + z 4 + z 5 )T(-z 2 - z 5 ) . (3.8) 

An auxiliary analytic regularization, by means of X\ and x 2 , is introduced into the powers 
of the propagators with the momenta p — ki and q — k 4 . The (lc-lc-lc) contribution can be 
obtained from (3.8) by the permutation x\ <-> x 2 . Each of the two (c-c-c) contributions is 
singular at x±,x 2 — > 0. The singularities are however cancelled in the sum. It is reasonable 
to start by revealing this singularity. One can observe that it appears due to the product 

T(2 - e - x 1 + z l + z 2 + z 3 + z 4 + 2 5 )r(— 2 + e + rr 2 - zi - z 2 - z 3 - z 4 - z 5 ) , (3.9) 

where the sum of the arguments of these gamma functions is x 2 — x±. 

So, the starting point is to take minus the residue at z 5 = —2 + e + x 2 — z 1 —z 2 — z 3 — z 4 and 
shift the integration contour correspondingly. The value of the residue is then symmetrized 
by x\ <-> x 2 . This sum leads, in the limit xi,x 2 — > 0, to the following fourfold MB integral: 



4 v ' ' (-sy+^t 2 r(-2e) 

_j_ /- + *°° /- +lo ° ^ d2 r(i + zQr(-i - e - ^ - s 2 )r(-* 2 )r(-e - * 3 )r(-* 3 ) 
(2myJ_ too "J_ ioo Al J r(i-^)r(i-^)r(i-2 e + ^ + ^ 2 + ^ 3 ) 

r(i + gi + ,2 2 + * 3 )r(i + e + 2: 2 + * 3 )r(-i -e + z 2 - z 4 )r(-i - z x - z 4 )r(i + z 4 ) 

r(-2 - 4e - - z 4 )r(-l + 2e - Zi - z 4 )r(3 + Zi + z 4 ) 
xT(-l + 6-^1-22-^3- ^)r(2 + zi + z 4 )r(-2 - e - - z 4 )r(2 - e + ^ + z 3 + z 4 ) 

x [2 7 + L + V(-3e) + V(-2e) - ^(1 + 3e) - ^(1 + zi + z 2 + 23) + ^(-1 - 21 - z 4 ) 

-V>(-2 - 4e - 21 - z 4 ) + V(-l + 2e - z 1 - z 4 ) + -e + z 2 - z 4 ) 

+ e - zi - z 2 - z 3 - z 4 ) + ^(2 - e + zi + z 3 + z 4 )\ , (3.10) 

where L = ln(s/t). 



r(3)b,c—c—c,res 
4 



( S ,t) 



r(-3e)r(l + 3e) 
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In the integral over the shifted contour in z 5 , one can set X\ = x 2 = to obtain the 
following fivefold integral: 

2e 3e7 r(l + 3e)r(-3e) 



T (3)b,c—c—c,int/ ,\ 

i\ (s, t) 



<- s y+3e(_ t y r(_ 2e ) 



X 



1 



+ioo 



r(l + ^)r(l + e + z 2 + Z3 )T(-z 2 )T(-e - z 3 )T(-z 3 ) 

r(i - z 2 )v(i - z 3 ) 



X 



-too j=1 

r(-i -e-z!- z 2 )T(i + Zl + z 2 + * 3 )r(-i - z x - z 4 ) 

T(l - 2e + Zl + z 2 + z 3 )F(-2 - 4e - z x - z 4 ) 
r(i + z 4 )r(e - z 2 - z 3 - z 5 )r*(-2 + e - - z 2 - z 3 - z A - z 5 ) 



r(l + e - ,2 2 - -23 - z 5 )T(l + e + z 2 + z 3 + z 5 ) 
xr(i + z 5 )r(-i - e + z 2 - z 4 )r(-2e + z 2 + z 3 + z 5 ) 

xT(2 - e + zi + z 2 + z 3 + z 4 + z 5 )T(-z 2 - z 5 ) , 



(3.11) 



where the asterisk on one of the gamma functions implies that the first pole is considered 
to be of the opposite nature. 

The evaluation of (3.10) and (3.11), in an expansion in e, is then performed according to 
the strategy characterized above. After the resolution of the singularities in e one obtains 
60 contributions where an expansion of the integrand in e becomes possible. Eventually, one 
reproduces the following leading asymptotic behavior: 

1 
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(3.12) 



3265920 

To compare eq. (3.12) with the complete result (B10), we use transformation formulae 
such as eq. (A7) to invert the arguments of the harmonic polylogarithms. The resulting 
quantities H aia2 „ Mn (l/x) with a n — 1 vanish as x — > oo. Logarithms are generated by 
the transformation; these logarithms combine with the ones already manifest in eq. (B10), 
yielding an expression in complete agreement with eq. (3.12). 
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IV. ITERATIVE STRUCTURE OF AMPLITUDES 

The iterative structure of the four-point MSYM amplitude found at two loops is [1, 13], 

Mf\e) = \{M^\e)) 2 + f^(e) Mj 1} (2e) + + 0(e) , (4.1) 

where 

f {2) (e) = -(( 2 + ( 3 6 + ( 4 e 2 + ■■■), (4.2) 
and the constant is given by 

C^ = -\d. (4.3) 

This relation can be verified by inserting the expansion (B5) for the planar double-box 

(2) (2) 

integral I\ into eq. (2.4) for M\ (e), and the expansion (B2) for the one- loop box integral 
1^ into eq. (2.3) for M^\e). Up through the finite terms in e, only harmonic polylogarithms 
with weights up to four are encountered (see appendix A). These functions can all be written 
in terms of ordinary polylogarithms if desired. Non-trivial cancellations between weight-four 
polylogarithms are needed to obtain eq. (4.1), strongly suggesting that the relation is not 
accidental and leading to the conjecture that an iterative structure exists in the amplitudes 
to all loop orders. 

As mentioned in the introduction, the relationship (4.1) is valid only through O(e ), 
i.e. near four dimensions, where MSYM is conformal and the AdS/CFT correspondence 
should be applicable. At 0(e 1 ), the difference between the left- and right-hand sides is an 
unenlightening combination of weight-five harmonic polylogarithms, not a simple constant. 

In order to search for a relation similar to eq. (4.1) at three loops, we have substituted the 
values of the integrals and I^ b , given in eqs. (B7) and (B9) respectively, into eq. (2.5) 
for M < f\e). We have also used the e-expansions of the one- and two- loop amplitudes through 
C(e 4 ) and 0(e 2 ) respectively, two further orders than required for the two-loop relation (4.1). 
(We cannot use eq. (4.1) to replace M^f 1 with because that equation is valid only 

through O(e ).) Thus we have obtained a representation of the amplitudes in terms of 
harmonic polylogarithms [24, 25] with weights up to six. Because harmonic polylogarithms 
with arguments equal to —t/s and —s/t both appear, we need to employ identities which 
invert the argument, of the type outlined in appendix A. 

Motivated also by the structure of the three- loop IR divergences described in ref. [26], we 
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have found the following iterative relation for the three-loop four-point amplitude, 



Mf(e) = -I[Mi 1) (e)]VMi 1) (e)Mi 2) (e) + /( 3 )(e)Mi 1) (3e) + C'( 3 ) + 0( e ), (4.4) 



where 



f (3 \e) = y C4 + e(6C 5 + 5C 2 Cs) + e 2 ( Cl C 6 + c 2 Q) , 



(4.5) 



and the constant is given by 



C (3) 



5 V / 17 2 , 



(4.6) 



'341 2 
-216 + 9 ( 

The constants c\ and c 2 are expected to be rational numbers. They do not contribute to 
the right-hand side of eq. (4.4) because of a cancellation between the last two terms, so they 
cannot be determined by our four-point computation. The reason we introduce them is to 
handle the subsequent generalization to the n-point MHV amplitudes. 



A. An ansatz for planar MHV amplitudes to all loop orders 

The resummation and exponentiation of IR singularities described by Magnea and Ster- 
man [30], and the connection to n-point amplitudes discussed by Sterman and Tejeda- 
Yeomans [26] (both of which we shall review shortly), together with the two- and three- loop 
iteration formulae, motivate us to propose a compact exponentiation of the planar MHV 
n-point amplitudes in MSYM at L loops. We propose that 



OO r OO 

M n = 1 + £a L M< L >(e) = exp £a l (fW(e)MW(U) + C« + E«\e)) 

L=l W =1 



(4.7) 



In this expression, the factor, 

0=^(4*0% (4.8) 

keeps track of the loop order of perturbation theory, and coincides with the prefactor in 
brackets in eq. (2.1). The quantity Mn\le) is the all-orders- in- e one-loop amplitude, with 
the tree amplitude scaled out according to eq. (2.2), and with the substitution e — > le 
performed. Each /^(e) is a three-term series in e, beginning at O(e ), 

f {l) (e)=ti l) + efi l) + e 2 fM. (4.9) 

The constants and are independent of the number of legs n. They are polynomials in 
the Riemann values ( m with rational coefficients, and a uniform degree of transcendentality, 
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which is equal to 21 - 2 + k for f k l \ and 21 for C {1) . The f [ k l) and C (l) are to be determined 
by matching to explicit computations. The En\e) are non-iterating O(e) contributions to 
the /-loop amplitudes, which vanish as e — > 0, En\o) = 0. 

Let us first see how eq. (4.7) is consistent with the results up to three loops discussed 
earlier in this section, by matching the left- and right-hand sides of the equations order-by- 
order in a. The one-loop case is very simple, since we only have to expand the right-hand 
side of eq. (4.7) to 0(a). It agrees with the left-hand side provided that 



c« = o, 



.10) 



That is, by definition we have absorbed the all-orders-in-e one-loop amplitude into M^\e). 
(It is possible that for n > 4 a nonzero value of En\e) could be more natural, given what 
is known about the structure of the one-loop amplitudes at 0(e) [56].) 

Next we expand eq. (4.7) to two loops, or 0(a 2 ). Using the one-loop result (4.10) to 
rewrite the 0(a) term in the exponential on the right-hand side of eq. (4.7) as Mn\e), we 
find that 

(6) = \ [M«(e)] 2 + f^\e) M«(2e) + + E^(e), (4.11) 

which is just the generalization of eq. (4.1) to n external legs. Evidence based on collinear 
limits in favor of this n-leg version, which we shall review in section V, was presented in 
ref. [1]; the values of /^(e) and given in eqs. (4.2) and (4.3) are independent of n. 

At the three- loop level, we also use the two-loop result (4.11) to rewrite the 0(a 2 ) term 
in the exponential on the right-hand side of eq. (4.7) as Mn\e) — \[Mn\t)} 2 . Matching 
both sides at 0(a 3 ) gives, 

M«(6)] 3 + M«( C ) {M^(e) - \ [M«( e )] 2 } 
{e)M^(?>e) + C^+E^\e) 

(e)l 3 + M«(c) Mf(e) + /< 3 >(c) M«(3e) + + E®(e) . (4.12) 



1 

6 

+ / 0), 
1 



For n — 4, this equation is equivalent to eq. (4.4), with the identifications (4.5) and (4.6) 
for /( 3 )(e) and . 

Eqs. (4.11) and (4.12) are special cases, for L = 2 and 3, of a more general L-loop iteration 
formula implied by eq. (4.7), 



M (L) = X W[M»(c)] + f {L \e)M^\Le) + + Ef\e) 



(4.13) 
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The quantities X^ = X^^Mn 1 ] only depend on the lower-loop amplitudes Mn'(e) with 
/ < L. For L = 2, 3, the values of X^ are, from eqs. (4.11) and (4.12), 

Xi 2) [M«Ke)} = \[M^]\ (4.14) 
XL 3) [M^(e)] = -\[m^]\m^M^. (4.15) 



(0, 



Now we establish eq. (4.13) for arbitrary values of L, and provide a convenient way to 
compute the functional Xn L \M^]. Taking eq. (4.13) as a definition of xi L \ we see that the 
full amplitude M n in eq. (4.7) can also be written as, 



M n = 1 + ^ = eX P E a "( M n L) - X( n ] ) 



-L=l 



(4.16) 



We need to show that the X^ only depend on the lower-loop amplitudes M„' with I < L. 
This result can be established inductively on L by comparing the 0(a L ) terms in the two 
Taylor expansions. The coefficient of a L on the left-hand side of eq. (4.16) is On 
the right-hand side, occurs explicitly in the a L term, and with the right coefficient to 

match the left-hand side. Every other term on the right-hand side depends only on M$ 
with I < L (using induction for those X- a L -* with V < L). But X^ must cancel all these 
other terms for the two Taylor expansions to agree; hence it also depends only on M$ with 
I < L. 

To solve eq. (4.16) for Xn L \ we take the logarithm of both sides, and look at the L th term 
in the Taylor expansion of that expression. We obtain, 



i=i 



(4.17) 



a L term 



Eqs. (4.13) and (4.17) are key equations; together they provide an explicit recipe for writing 
the L-loop amplitude in terms of lower-loop amplitudes, plus constant remainders. 

From eq. (4.17) we can easily recover eqs. (4.14) and (4.15), as well as obtain, for example, 
the next two values of X^\ 

X^[M«\e)\ = \[m^\- [m«] 2 Mf+M«Mf + \[m^}\ (4.18) 



(0/ 
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Mi 3 ) 



(4.19) 
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Note from eq. (4.7) that / w (e) appears multiplied by M n (/e) , which has poles beginning 
only at order 1/e 2 . Hence we can absorb any 0(e 3 ) and higher terms in /^(e) into the defi- 
nition of the non-iterating contributions E n l \e). However, the 0(e 2 ) terms in /^(e), namely 
ff \ cannot be removed because is asserted to be independent of n. This statement can 
only be true for one choice of ffp; shifting that value induces a shift proportional to n in 
C®, because M^\h) oc n/e 2 . The value of / 2 (0 can be determined by computing an /-loop 
amplitude with n > 4, or else the /-loop splitting amplitude (which may be simpler), as 
reviewed in section V. 



B. Infrared consistency of ansatz 

In this subsection we discuss the consistency of the exponentiated L-loop ansatz (4.7) 
with the resummation and exponentiation of IR divergences [29], following the analysis of 
Magnea and Sterman [30], and of Sterman and Tejeda-Yeomans [26]. 

A general n-point scattering amplitude can be factorized into the following form, 

M„ = J (Q^, a s (n), ejx5 a s (/i), e^j x h n a s (fx), e^j , (4.20) 

where J is a jet function, S a soft function, and h n a hard remainder function which is finite 
as e — > 0. Also, fx is the renormalization scale, and Q some physical scale associated with 
the scattering process for external momenta k{. 

Both M. n and h n are vectors in a space of possible color structures for the process, 
and S is a matrix. However, we shall work in the leading-color (planar) approximation, in 
which there is no mixing between the different (color-ordered) color structures. Hence S is 
proportional to the identity matrix. As pointed out in ref. [26], S is only defined up to a 
multiple of the identity matrix, so we can absorb it into the jet function J at leading color. 
Figure 5 illustrates that, at leading color, soft exchanges are confined to wedges between 
color-adjacent external lines, for example the lines i and i + 1 in the figure. We also consider 
adjoint external states, such as gluons. By examining the case n — 2, it can be seen that the 
wedge that is being removed in the figure represents half of the IR singularities of a Sudakov 
form factor [29]; that is, a color-singlet object (such as a Higgs boson) decaying into two 
gluons. We denote this matrix element as M^ 9 ^(s iii+ i/ /i 2 , a s (/x), e). 

Because MSYM is conformally invariant (the f3 function vanishes), a s may be set to a 
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FIG. 5: Infrared structure of leading-color scattering amplitudes for particles in the adjoint rep- 
resentation. The straight lines represent hard external states, while the curly lines carry soft or 
collinear virtual momenta. At leading color, soft exchanges are confined to wedges between the 
hard lines. 

constant everywhere. Thus the leading-color IR structure of n-point amplitudes in MSYM 
may be rewritten as, 



M n = n 



1=1 



M 



[9S-1] [ 



2 i ®-si £ 



1/2 



(4.21) 



where h n is no longer a color-space vector. 

For a general theory, the Sudakov form factor at scale Q 2 can be written as [30] 



M [ "~ 



exp 



2 

l6 a \a s {n), e) + <?M (-1, a s (^, a.(ji), e) , e) 

(4.22) 



1 [» dfi 2 [g] /_ /fx 2 



where 7^ denotes the soft or (Wilson line) cusp anomalous dimension, which will produce 
a 1/e 2 pole after integration. The function K,w is a series of counterterms (pure poles in e), 
while includes non-singular dependence on e before integration, and produces a 1/e pole 
after integration. 

In MSYM, a s (ii) is a constant, and the running coupling a s (/x 2 //i 2 , a s , e) in 4 — 2e di- 
mensions has only trivial (engineering) dependence on the scale, 



This simple dependence makes it very easy to perform the integrals over £ and fl. 



(4.23) 
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Following refs. [26, 30], we expand K} 9 \ 7^, and in powers of a s , 

00 1 

6^(-l,a,(^, Q „ £ ), £ ) =E«'@>, 



(4.24) 
(4.25) 
(4.26) 



where a is defined in eq. (4.8) and the hats are a reminder that the leading-iV c dependence 
has also been removed in eqs. (4.24), (4.25) and (4.26). That is, the perturbative coefficients 
(defined with expansion parameter a s /(2n)) have a leading-color dependence on N c of, 



(4.27) 



We can suppress the [g] label because the N = 4 MHV amplitudes are all related by 
supersymmetry Ward identities [57], so that the corresponding functions for external gluinos, 
etc., are the same as for gluons. Equation (4.24) follows from solving eqs. (2.12) and (2.13) 
of ref. [30] in the conformal case ((3 = 0). In this case, K)& contains only single poles in e, 
which are simply related to 7^. 
The integral over Q is very simple, 



d_e 

Jo t 



i 2 ^ le \-Q 2 J ' 



1=1 



(0 




(4.28) 



The first integral over 7^ gives, 



-(0 
Tk ■ 



(4.29) 



Adding the K$ term to 1/2 of eq. (4.29), using eq. (4.24), we see that the "—1" is 
cancelled. Then the integral over £ is properly regulated, and evaluates to 



OO ; 



Combining this result with eq. (4.28) gives 



" 2 V< 7<? 



(4.30) 



1=1 



(4.31) 
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We need eq. (4.31) for a neighboring pair of legs i,i + 1 in the n-point amplitude, so 
that Q 2 should be replaced by the invariant s iti+ i. Taking the product over all i, eq. (4.21) 
becomes 



M n = exp 



i=i 



X h r , 



We may rearrange this a bit, to give 



M n = exp 



i=i 



x h n , 



where /^(e) is defined in eq. (4.9), with the identifications, 



Jo 



(0 



(0 



1 .(o 
l -G il) 



(4.32) 



(4.33) 



(4.34) 
(4.35) 



and 



-i -i 9 

1 1 ^^^/ /X 



2 ^ ^II V ~ S M+1 



(4.36) 



Here fo n differs from /i n by a finite shift, due to the 0(e 2 ) terms in /^(e), which we introduce 
to help make contact with the exponentiated ansatz (4.7). Using /^(e) = 1 and eqs. (4.2) 
and (4.5), we may read off the first few loop orders, 



7? =4, 
# = 22 C 4 , 



(4.37) 



and 



= o, 



C(2) 



-Cs, 

4Cs + y C2C3 



(4.38) 



The quantity In\e) is a function that captures the divergences of the planar one-loop 
n-point amplitudes in MSYM [27], after extracting the leading- N c dependence as in ref. [13]. 
The W defined in eq. (8.9) of ref. [13] contained a prefactor of e e ^( 1 )/r(l — e), following 
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conventions of Catani [28]. Here we adopt a convention closer to that of Sterman and Tejeda- 
Yeomans [26], without such a prefactor. The difference between eq. (8.9) of ref. [13] and 
eq. (4.36) above is a finite quantity, because e~ € ^ /T(l — e) = 1 + 0(e 2 ). Finite remainders 
will differ between the two conventions. 

Starting from eq. (4.33), and using the fact that the difference between Mn\le) and 
In\le) is finite, we can reshuffle the finite terms once more to obtain, 



M n = exp 



• oo 

£«< (f {l) (e)MM(le)+h«\e)) 



(4.39) 



■i=i 

We have moved the hard function into the exponent without loss of generality, because we 
allow for a new function h$ at each order I. 

Finally we compare the exponentiated ansatz (4.7) with the exponentiation of the IR 
divergences (4.39). We see that they agree if we identify 

h$>(k i ,e) = CW + Eg>(e). (4.40) 

In some sense, the content of the iterative structure of planar MSYM, beyond the level of 
consistency with IR resummation, is that the (suitably-defined) exponentiated hard remain- 
ders hn\ki,e) approach constants, independent of the kinematics and of n, as e — > 0, since 
E%\e) is of 0(e). 



C. Finite Remainders 

Next we shall obtain iterative formulae for two series of functions: the 1^ governing IR 
divergences for the L-loop n-point planar amplitudes, and the representing the finite 
remainders of the amplitudes, after subtracting these divergences. The formulae will be very 
similar in form to the full amplitude relation (4.13). 

Following the structure uncovered explicitly at one, two, and three loops [26, 27, 28], we 
define the finite remainder for the L-loop amplitude by writing 

Mi L He) = J2li L - l \e) AfW(e) + *f >(e) , (4.41) 

1=0 

or ^ 

i* L) (e) = MW-£#-0M«, (4.42) 

1=0 
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where Mn^ = 1. We insert the iteration formula (4.13) for the first term, M-i L \ on the 
right-hand side of eq. (4.42). We split M^(Le) -> l£\Le) + F^\Le) in this formula. For 
the lower-loop amplitudes, Mn\ we recursively substitute in the finite-remainder formula 
for smaller values of /, 

Mg>= f^it k) M<P+Fg>. (4.43) 



fc=0 



At this point, the expression for F^\e) is a polynomial in ijp and Frp , which has the special 
property that there are no mixed I-F terms. (If there had been such terms, it would have 
signaled an inconsistency.) We can remove the pure-/ terms by choosing 1^ to cancel them. 
The resulting finite expression gives Fn L \e) as a polynomial in the lower-loop Fn\e). 

We find that the solutions for li L \e) and Fn L \e) are expressible in terms of the same 
function defined in eq. (4.17), but where the role of M n is played instead by — I n and 
F n , respectively: 



= -XL L \-m}+f {L) (e)lU(Le), 
Fi L) (e) = Xi L \F«\e)] + f^\e)F£\Le) + C^+Ei L \e) 



(4.44) 
(4.45) 



The Taylor expansion (4.17) can be used to evaluate eqs. (4.44) and (4.45) to any desired 
loop order. 

Because the form of eq. (4.45) for F^\e) is completely analogous to the iteration for- 
mula (4.13) for the full amplitude M^ L \e), we see that the finite remainders can be expo- 
nentiated as, 



T n {e) = l + Xyi^e) =exp £ a 1 (/»(c)FW(Zc) + C« + S«(c)) 

L=l H=l 

Letting e — > 0, we have 

OO r OO 

F n (0) = 1 + ^V^f >(0) = exp £ a 1 (/<'>*?> (0) + C«) 

L=l >-l=l 

where 

/o ( °-/ ( °(0). 



(4.46) 



Using eq. (4.34) we may then rewrite this as, 



.F„(0) = exp 



\ikFW(0) + C 



(4.47) 



(4.48) 



(4.49) 
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The soft anomalous dimension is 

oo 

Ik = J2^K al = 4a-4C 2 a 2 + 22C 4 a 3 + --- , (4.50) 
1=1 

where we used eq. (4.37). Similarly, from eqs. (4.3), (4.6) and (4.10), we have, 

oo 1 

C = J2C (l) a l = --C 2 2 a^ + 
i=i 

As mentioned below eq. (4.6), the rational numbers C\ and c 2 are yet to be determined. The 
resummation (4.49) of the finite remainders of the MHV amplitudes, as a consequence of 
the exponentiated ansatz (4.7), is one of the key results of this paper. 

For the Fn\o), the argument le in eq. (4.45) has disappeared as e — > 0. Hence we can 
recursively substitute back to obtain formulas solely in terms of Fn\o). Equivalently, we 
can series expand the exponential in eq. (4.47) or (4.49). The first few values are, 

Fi 2 \0) = ^(O)) 2 + /o (2) Fi 1} (0) + C (2) , (4.52) 
FfW = -l(Fi 1 \0)y + FW(0)F^(0) + fl ) 3) FW(0) + C^ (4.53) 
= ^(V ) (0)) 3 + /f(Fi 1) (0)) 2 + (/W +C«)^(0) + CW, (4.54) 

^i 4) (o) = ^Ho)) 4 -(^ 

+/o (4) ^ 1} (0) + ^ (4) (4.55) 
= ^(^(0)) 4 + ^o (2) (^ 1} (0)) 3 + + 2/ (3) + C^) (i^(0)) 2 

+ (/o (4) + ti 2) C {2) + C^)FW(0) + CW + l -[C^f . (4.56) 

Thus, starting from the ansatz (4.7), we have succeeded in expressing the n-point L-loop 
MHV finite remainders directly in terms of the one-loop finite remainders. 

We remark that the two-loop result (4.52) differs slightly (in the constant term) from 
the corresponding eq. (16) for n = 4 in ref. [1]. The reason is that the definition of the 

"(2) 

two- loop divergence used there, interpreted in terms of In , does not obey eq. (4.44) for 
L = 2, but differs from that I n 2 ^ by a finite (O(e )) amount. The definition we use here is 
more convenient because of its simple generalization to higher loops. 

The one-loop finite remainders F^(0) for the MHV amplitudes in MSYM were evaluated 
for all n in ref. [9], using the unitarity method. Modifying those results to the conventions 



/341 2 \ / 17 2 \ 
V216 + 9 Cl J C6+ (-Y + 9 C2 J C: 



a 3 + 



(4.51) 
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of this paper, the finite terms are explicitly, for all n > 5, 

F£\V) = \j^9n^ (4.57) 

where 

L"/2J-1 , _.[r] \ / _.[r] v o 



and where \_x\ is the greatest integer less than or equal to x. Here tf = (fcj H h /c i+r _i) 2 

are the momentum invariants, so that ij 1 ' = and tf^ = s^j+i. (All indices are understood 
to be modn.) The form of D ni and L n>i depends upon whether n is odd or even. For 
n = 2m + 1, 

m-l / jMjH- 2 ! 



r=2 ^ J i ' 

1 / _/[ m l \ / -t [m] \ 

= -2 In (-pi— In -gl , (4.60) 

whereas for n = 2m, 

m-2 , JrUr+2] , , [m-l],[m+lK 

jD2mt = _y Li2 (i_ ^-i )_i Li2 (i_ ^ ), (4.6i) 

1 / _ / M \ / _ f H \ 



4 V _+H / V _+H / ' 
For n = 4 the above formula does not hold, but the finite remainder is simply, 

if } (0) = ihi 2 (3+ 4 C2- (4-63) 

Assuming that the exponentiated ansatz (4.7) holds, then the exponentiated finite re- 
mainders .F n (0) given in eq. (4.49) are completely determined to all loop orders, in terms of 
the one-loop remainders Fn\o) just presented, plus the series of constants 7^ and C given 
in eqs. (4.50) and (4.51). 

V. COLLINEAR BEHAVIOR AND CONSISTENCY OF ALL-ra ANSATZ 



In this section we discuss the consistency of the n-point iterative ansatz (4.13) with the be- 
havior of amplitudes under factorization. In a supersymmetric theory, MHV amplitudes have 
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no multi-particle poles; the residues vanish by a supersymmetry Ward identity [57]. This 
property is manifest in our ansatz, because neither the tree amplitude An^' MHV (1, 2, . . . , n) 
nor the one-loop amplitude M^\{le) contain such poles. Hence only factorizations as pairs 
of momenta become collinear need to be analyzed. 

In general, color-ordered amplitudes Ai L ^(l, 2, . . . , n) satisfy simple properties as the mo- 
menta of two color- adjacent legs k h k i+1 become collinear, [9, 12, 41, 58, 59], 

L 

Af\. . . (i + 1)**V ••)—>££ Split® (z;i\ (i + l) x ^)Ai L _- l \. . . , P*. . .) . (5.1) 

The index / sums over the different loop orders of contributing splitting amplitudes Split 
while A sums over the helicities of the intermediate leg kp = —{ki + fcj+i), and z is the 
longitudinal momentum fraction of fcj, ki ~ —zkp. The splitting amplitudes are universal 
and gauge invariant. 

The tree-level splitting amplitudes Split ^ are the same in MSYM as in QCD. At loop 
level, the MSYM splitting amplitudes are all proportional to the tree-level ones. The pro- 
portionality factors depend only on z and e, not on the helicity configuration, nor (except for 
a trivial dimensional factor) on kinematic invariants [9]. It is thus convenient to write the 
L-loop planar splitting amplitudes in terms of "renormalization" factors r^\e; z, s), defined 
by 

Split5! p (l Al ,2 A2 ) = 4 L) (6;z, S )SplitL°i p (l Al ,2 A2 ), (5.2) 
where s — (ki + k 2 ) 2 . 

Using eqs. (5.1) and (5.2), we see that the amplitude ratios Mn L \e) = /A^ behave 
in collinear limits as, 

M n\ e ) - M^+r^ie), (5.3) 
Mj?He) - M^e) + rf (ejM^e) + rf(e) , (5.4) 

or at three loops, 

Mf(e) - M^(e) + ^(e) M^e) + rf^M^e) + rf{e) , (5.5) 
where rf\e) = 1 and we have suppressed all functional arguments except for e. 
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The one-loop splitting amplitudes in MSYM have been calculated to all orders in e, with 
the result [58], 

,2 



(e;z,s) 



c r ( \l 



e 2 v— s 



Tie 



sin [Tie) 



l-z 



z 



2fc+l 



Li 



2fc+l 



-Z 



fc=0 



l-z 



(5.6) 



where Li„ is the n-th polylogarithm (defined in eq. (A5)), and 

e^ r(l + e)r 2 (l-e) 

Cr = T r(i - 2e ) • (5 ' 7) 

In refs. [1, 13], the two-loop splitting amplitudes in MSYM were computed through G(e ) 
using the unitarity method as described in ref . [59] . The result of this computation is a very 
simple formula, expressing the two-loop splitting amplitude in terms of the one-loop one, 

4 2) (e; z, s) = \ (rf{e- z, s)) 2 + f^\e) rf(2e; z, s) + 0{e) , (5.8) 

where f^(e) is given in eq. (4.2). (This result was actually obtained before the iterative 
relation (4.1), and motivated its discovery.) 

The consistency of the n-point ansatz (4.13) for L = 2 (with X„ given by eq. (4.14)) 
may be easily confirmed using these splitting functions [1]. Inserting the collinear behavior 
of the one-loop amplitudes (5.3) into the right-hand side, we obtain, 

= M^{e) + rf{e)M^_ 1 {e) + rf{e) , (5.9) 

where we used eqs. (4.13) and (5.8) to rearrange the expression into the required form (5.4) 
for correct two-loop collinear behavior. Since there are no multi-particle poles in the MHV 
case, eq. (5.9) confirms that the ansatz (4.13) has the correct factorization properties in all 
channels at two loops. 

Similarly, we can require that the ansatz (4.13) is consistent with collinear factorization 
beyond two loops, and thereby obtain an iterative ansatz for the planar L-loop splitting 
amplitudes in MSYM. For L = 3, the ansatz reads, 

Mf( e ) = -i[MW(e)] 3 + AfW( e )MW (e) + f^\e) (3 e) + + 0(e) . (5.10) 

If we insert the properties of one- and two-loop amplitudes (5.3) and (5.4) into the collinear 
limit of the right-hand side of eq. (5.10), we obtain, 



Mf(< 



M« 1 (6)+4 1) (e) 
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+ f i3) (e) (m^c) + r«(3c)) + C< 3 > + 0(e) . (5.11) 

After rearranging terms, we can get consistency with eq. (5.5), provided that the three- loop 
splitting function obeys, 

rf(e) = -l(rf\e)) 3 + r< f 1) (e)r<? ) (e) + f< 3 \e)r^(3e) + (9(e) . (5.12) 

By repeating this exercise at L loops, and collecting the terms that are independent of 
^l-i) we see that the relation, 

rf(e) = X^[4\e)} + /W(e) rf(Le) + 0{e) , (5.13) 

is the one required for consistency. (We have dropped the subscript n from because it 
is out of place here, but it is the same function of lower- loop quantities defined in eq. (4.17).) 
In other words, the L-loop splitting amplitude functions rif obey exactly the same type 
of iterative relation as the scattering amplitudes M^ L \ but without the "inhomogeneous" 
constant terms C^ L \ Because the one- loop splitting amplitude r^(Le) begins at order 
0(e~ 2 ), the relation (5.13) allows the 0(e 2 ) coefficient of f <yL \e), namely f^ L \ to be extracted 
from the 0(e°) term in the L-loop splitting amplitude. 



VI. ANOMALOUS DIMENSIONS AND SUDAKOV FORM FACTORS 

The soft anomalous dimension ^ K controlling the 1/e 2 IR singularities of the loop am- 
plitudes arises from an edge of phase space, the Sudakov region, where a hard line can 
only emit soft gluons. In the loop amplitudes these gluons are virtual, of course, but they 
are related to real soft-gluon emission by the cancellation of infrared poles in infrared-safe 
cross sections for Sudakov-type processes [29, 31]. For example, the splitting kernel Pu(x) 
describes the probability for a parton i to split collinearly into a parton of the same species 
i, plus anything else, where the second parton i retains a fraction x of the longitudinal 
momentum of the first parton i. In the limit x — > 1, this splitting kernel is dominated by 
soft-gluon emission, and has the form, 

P«(x) -> + B(a s )8(l -x) + ..., as a; -> 1, (6.1) 

[i - x) + 
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where A(a s ) is related [35] to the soft (cusp) anomalous dimension by, 



(6.2) 



The splitting kernel is related by a Mellin transform to the anomalous dimensions of leading- 
twist operators of spin j, 



Jo 

Thus the soft anomalous dimension also controls the large-spin behavior of these anomalous 
dimensions [60], 



where here we take 7 e as Euler's constant. 

KLOV [38] have made a very interesting observation: the anomalous dimensions of MSYM 
may be extracted directly from the corresponding anomalous dimensions of QCD [36], by 
keeping terms of highest "transcendentality" . Recall that for the case of the soft anomalous 
dimensions (large j limit), the transcendentality weight is simply n for ( n . (Although the 
QCD anomalous dimensions are computed in the MS regularization scheme, whereas for 
MSYM the DR [45] or FDH [44] schemes are needed to preserve supersymmetry, the scheme- 
dependent terms drop out because they are of lower transcendentality.) Although there is no 
proof of KLOV's prescription for extracting the MSYM anomalous dimensions from QCD, 
there are good reasons to believe that it is true [13, 38]. 

Here we provide further evidence for the prescription, by confirming the large-spin be- 
havior of the leading- twist anomalous dimension. We compare the KLOV result, given in 
eqs. (18)-(20) of ref. [38], against our evaluation of the same quantity from the IR diver- 
gences of the three-loop four-point amplitude. Note that their normalization convention 
for anomalous dimensions has an opposite overall sign from ours (which follows ref. [36]). 
Also taking into account factors of 2 from the different a s expansion parameter, and from 
eq. (6.4), we obtain from eqs. (18)-(20) of ref. [38], 




(6.3) 



70') = 2"7*K) 0*0) + 7e) " B(a s ) + 0(\n(j)/j) , 



(6.4) 




AN, 




-<2 N 2 C 



c ' 



(6.5) 




22C 4 iV c 3 , 



which agrees perfectly with our results (4.27) and (4.37). 
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We remark that the strong-coupling, large- N c limit of the soft anomalous dimension 7^ 
has been obtained, using the AdS/CFT correspondence and classical supergravity meth- 
ods [61]. An approximate formula interpolating between the weak and strong-coupling 
limits has also been constructed [38, 62]. 

The coefficient Qq\ which controls the 1/e singularity, may be extracted [32] from a 
fixed-order computation of the form factor at / loops. For example, the two-loop quark form 
factor in QCD was computed in ref. [63]. From eqs. (21)-(22) of ref. [26], if we follow the 
KLOV procedure and keep the maximal transcendentality terms (£3 at two loops) in order 
to convert the QCD results into MSYM results, we have, 

si" = 0, 

6f = -6 A? ■ (6.6) 

(2) 

We have multiplied Qq in eq. (22) of ref. [26] by a factor of 4 to account for the different 
normalization conventions used here. These results agree with our eqs. (4.27) and (4.38). 
Although the QCD form factors have not yet been computed at three loops, we may use 
our results, together with the observation of KLOV, to predict the leading-transcendentality 
contributions for QCD, 

^ 3) = (4C 5 + yC 2 C 3 )iV c 3 , (6.7) 

after the group theory Casimirs have been set to the values C F = C A = N c . (At three 
loops, no other Casimirs can appear, so there are no subleading-color corrections to this 
leading-transcendentality prediction.) 

VII. CONCLUSIONS AND OUTLOOK 

In this paper we have provided strong evidence supporting the conjecture [1] that the 
planar contributions to the scattering amplitudes of MSYM possess an iterative structure. 
This result is in line with the growing body of evidence that gauge theory amplitudes in 
general, and those of MSYM in particular, have a much simpler structure than had been 
anticipated. 

Our evidence of iteration is based on a direct evaluation of the planar three- loop four-point 
amplitude of MSYM. The loop integrands for this amplitude were obtained [7, 8] using the 
unitarity method [9, 10, 11, 12, 13]. This method ensures that simple structures uncovered at 
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lower loop orders (including tree level) in turn feed into higher loops. (It also underlies much 
of the recent progress at one loop [19].) In order to evaluate the required three-loop integrals, 
we made use of important recent advances in multi-loop integration [14, 21, 22, 23]. The 
integrals are expressed in terms of well-studied harmonic polylogarithms [24, 25], making it 
straightforward to confirm the three-loop iteration. A rather intricate set of cancellations is 
required, amongst the harmonic polylogarithms, and between different loop integral types 
contributing to the amplitudes. 

Using our explicitly computed four-point amplitudes as a springboard, the known struc- 
ture of infrared singularities to all loop orders [26, 30], and the required factorization prop- 
erties of amplitudes, we constructed the ansatz for the resummed n-point all-loop MHV 
amplitudes given in eq. (4.7). After subtracting the IR divergences, the all-loop finite re- 
mainders (4.49) are given in terms of known one-loop n-point finite remainders, as well as two 
coefficients, one of which is the large-spin limit of the leading-twist anomalous dimensions. 

Very interestingly, the same set of leading-twist anomalous dimensions has recently been 
linked to integrability of MSYM by Beisert, Kristjansen and Staudacher [34, 39]. With 
the assumption of integrability, Staudacher [39] has reproduced the leading-twist anomalous 
dimensions at three loops for spin j up to 8. These anomalous dimensions were previously 
obtained by Kotikov, Lipatov, Onishchenko and Velizhanin [38] from the QCD results of 
Moch, Vermaseren and Vogt [36]. (Quite recently, Staudacher's Bethe ansatz analysis has 
been extended to extremely high spins, the region relevant here, confirming the prediction of 
KLOV for even values of j up to 70 [64].) If one were able to push this method to higher loop 
orders, and arbitrarily large spins, it would give very directly the soft anomalous dimensions 
appearing in our all-loop exponentiation of the MHV scattering amplitudes. 

Besides confirming the iterative structure of the scattering amplitudes, our paper provides 
non-trivial confirmation of the form of the three-loop divergences predicted by Sterman 
and Tejeda-Yeomans [26]. It also provides supporting evidence for a number of ansatze 
appearing in a variety of papers. In particular, we confirm, in the high spin limit, the inspired 
ansatze of KLOV, and (via KLOV) of Beisert, Kristjansen and Staudacher, for obtaining 
the leading-twist anomalous dimensions in MSYM. By making use of KLOV's link to QCD, 
via the degree of transcendentality, our work also checks indirectly a small piece of the 
three-loop splitting kernels in QCD, or equivalently the anomalous dimensions of leading- 
twist operators, computed by Moch, Vermaseren and Vogt [36]. The integrand [7, 8] used 
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in the computation of the planar three-loop four-point amplitude has not been completely 
proven, but the match between its IR singularities and the formulae of Sterman and Tejeda- 
Yeomans, plus the demonstration of its iterative structure through the finite terms as e — > 0, 
leaves little doubt as to its veracity. 

The properties found here and in ref. [1] bring up the possibility that the entire perturba- 
tive series of planar MSYM is tractable. The apparent simple structure of the MHV all-loop 
amplitudes suggests that a loop-level twistor string interpretation will be found [16, 65]. 
It would be important to first identify the precise symmetry responsible for this structure. 
A more complete understanding of the iterative structure of the amplitudes should lead to 
important insights into quantum field theory and the AdS/CFT correspondence. 

Note added 

Since the first version of this paper came out, an interesting paper has appeared [66], 
containing a technique for computing large classes of terms for multi-loop MSYM amplitudes 
with many external legs, which may shed further light on the iterative relations discussed 
here. Also, the prediction (6.7) for the leading-transcendentality terms in Qq for QCD has 
now been confirmed [67]. 
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APPENDIX A: HARMONIC POLYLOGARITHMS 



We express the amplitudes in terms of harmonic polylogarithms [24], which are general- 
izations of ordinary polylogarithms [68] . Here we briefly summarize some salient properties. 
A more complete discussion is given in ref . [24] . Recipes for numerically evaluating harmonic 
polylogarithms may be found in ref. [25]. 

The weight n harmonic polylogarithms H aia2 .. Mn {x) = H(ai,a2, ■ ■ ■ ,a n ;x), with Oj G 
{1, 0, — 1}, are defined recursively by, 

H aia2 ...a n (x) = I dtf ai (t)H a2 ... an (t), (Al) 

JO 



where 



f±i(x) = T ^, fo(x) = -, (A2) 
1 =F x x 

H ±1 (x) = tH 1 , H (x)=\nx, (A3) 
and at least one of the indices is non-zero. For all Oj = 0, one has 

H 0fi _ (x) = -\n n x. (A4) 

71/. 

If a given harmonic polylogarithm involves only parameters a; L = and 1, and the number 
of these parameters (the weight) is less than or equal to four, it can be expressed [24] in 
terms of the standard polylogarithms 



T . . . ^ z j f z dt T . . , 

3=1 3 J ° 1 

Li 2 (z) = - f^ln(l-t), (A5) 
Jo 1 

with n — 2, 3, 4, and where z may take the values x, 1/(1 —x), or —xj (1 — x). (For n < 4, not 
all of these values are required, due to identities.) Here we need only G {0, 1}, but weights 
up to six. In the Euclidean region for the planar four-point process, namely s < 0, t < 0, 
u > 0, with the identification x = —t/s, the argument x of the harmonic polylogarithms 
will be negative. 

The harmonic polylogarithms are not all independent; they are related by sets of identi- 
ties [24]. One set of identities, derived using integration by parts, 

H ai a 2 ...a p o(x) = \nxH aia2 __, ap (x) - H 0aia2 __, ap (x) - H ai0a2 ...a p (x) H aia2 _„ a p (x) , (A6) 
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allows one to remove trailing zeroes from the string of parameters a^. The remaining 
H ai a 2 ...an( x ) with a n = 1 are well-behaved as x — > 0; in fact they all vanish there. 

Because the integrals appear in the MSYM amplitudes with arguments (s, t) and (t, s), we 
need a set of identities relating harmonic polylogarithms with argument x = — t/s to those 
with argument y = —s/t = 1/x. As explained in ref. [24] (see the discussion near eqs. (55) 
of that reference), we may construct the required set of identities by induction on the weight 
of the harmonic polylogarithms. For the first few weights, in the region — 1 > x > 0, and 
letting L = In (s/t) = ln(— 1/x), we have, for example, 

H^y) = H x (x)-L, 
HoM = -H 0A (x)-±L 2 -^, 

H 1A (y) = H lil (x)-H 1 (x)L + h 2 , 

7T 2 1 

H ,o,i(v) = H 0fi ,i(x) --L--L 3 , 

#o,i,i(y) = H 0A1 (x)-H 0A!l (x) + H 0A (x)L + ^L 3 + C3, 

H lfiil (y) = -2H 0:OA (x) + 2H 0:1A (x)-2H , 1 (x)L-^-H 1 (x)-H 1 (x)H 0:1 (x) 

11^ 1 

- -H^L 2 - -L 3 + yL + H 0A (x)L + -L 3 - 2 C 3 , 
#i,i,i(y) = H 1AA (x) - H ltl (x)L + l -H x [x)L 2 - h 3 . (A7) 



APPENDIX B: INTEGRALS APPEARING IN FOUR-POINT AMPLITUDES 

In this appendix we collect various integrals that are needed as well as their values in 
terms of harmonic polylogarithms. We quote the results in the Euclidean (w-channel) region, 
s, t < 0. The analytic continuation to other physical regions is discussed in refs. [24, 25]. 

1. One-loop integrals 



Consider the (conveniently normalized) one-loop scalar box integral, depicted in fig. 6, 

d d p 

p 2 (p - ki) 2 (p — ki — k 2 ) 2 (p + k^) 2 



ll%,t) = -ie^~ d ' 2 ( ^ ^ , . (Bl) 
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FIG. 6: The one-loop box integral. 



The value of this integral, with x = —t/s, L = ln(s/i), is 



Cj(x,L) 



with 

c 2 = 

Cl = 

c = 

C-l = 
C-2 = 



4, 
2L, 

7 1 34 

• 2 #i(:r) + 2if o ,o,i(a0 - t^ 2 L + 2// ,i(^)^ + H^L 2 - -L 3 - — C 3 , 

O 3 3 



n'H^x) + 2^0,0,1 W - R n + 2H o,i{ x L + HiWL* - -L" - —Q , 
-2# 1A0 ,i(z) - 2H 0fijlil (x) - 2H 0>lfitl (x) - 2H 0fifijl (x) - 2H 0ihl {x)L 

- 2H 1>0tl (x)L + H 0tl (x)L 2 - H 1A (x)L 2 + ^H^L 3 - 

1 _ . 9fl 



^ 3 - -L 4 

- 7r 2 iJ 1;1 (x) + tt 2 H\ (x}L - i^L 2 + 2^(x)C - yL( 3 - ^vr 4 , 
2if 1,0,0,0,1(2;) + 2^1,0,0,1,1(0:) + 2^1,0,1,0,1(0:) + 2^1,1,0,0,1(0:) + 2H 0fififijl (x) 

,0,0,1,1(3;) ,0,1,0,1(3;) + 2^0,0,1,1,1(3;) + 2^0,1,0,0,1(3;) + 2^0,1,0,1,1(2;, 

+ 2^0,1,1,0,1^) + 2H 0:1A:1 (x)L + 2H lfi ^(x)L + 2H lilfi , 1 (x)L - H^{x)L 2 

- H lfijl (x)L 2 + H lxl (x)L 2 + ±H 0il (x)L 3 - \h^{x)L 3 + ^i(x)L 4 - 1 

- ^H , 0>1 (x) + 7r 2 //i,i,i(x) - ^7T 2 // ,i(2;)L - n 2 H hl (x)L + ^k 2 H x {x)L 



C-3 



-— L 5 
20 



6 

5 

36' 



- 4vr 2 L 3 + ^tt 2 C 3 - 2//i,i(x)C 3 + 2E x {x)L^ - L 2 ( : 
10 



(B2) 



C-4 = 



7 4, 1 4„, , 134 

1,0,1,1(3;) — 2ifo,o,o,i,o,i(3;) — 2i?o,o,o,i,i,i(3;) — 2ifo,o,i,o,o,i(2;) 

{.,.} Oil (^\ _Off„. . (~\ _ o w_ . . ^ 



,0,0,0,0,1(3;) — 2iJ ,o,o,o 
~~ 2H 00 i : 
- 2# ,i,o, 



,0,0,0,0,1 (, 3; ; — ^-00,0,0,0,1,11,3;; — zn 0: o,o,i,o,i{ x ) — ^-"0,0,0,1,1,1 \ x ) ~ ^-"0,0,1,0,0,11,2;; 
0,0,1,0,1,1(2;) — 2iJ ,o, 1,1,0,1 ( x ) ~~ 2^0,0,1,1,1,1(2;) — 2iJ , 1,0,0,0, 1(2;) — 2^0,1,0,0,1,1(2;) 
0,1,0,1,0,1(2;) - 2^/0,1,0,1,1,1(2:) - 2^0,1,1,0,0,1^) - 2// ,i,i,o,i,i(2;) - 2^0,1,1,1,0,1 W 
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FIG. 7: The two-loop double-box integral. 

- 2H lj0fij0fijl (x) - 2^1,0,0,0,1,1 - 2 #i,o,o,i,o,i(z) - 2#i, ,o,i,i,i(z) - 2#i, ,i, ,o,i(z) 

- 2^1,0,1,0,1,1 (») - 2#i, ,i,i,o,i(z) - 2#i,i, ,o,o,i(z) - 2#i,i, ,o,i,i(z) - 2#i,i, ,i, ,i(z) 

- 2^1,1,1,0,0,1 (^) - 2# ,i,i,i,i(z)£ - 2^i, ,i,i,i(x)L - 2//i,i, ,i,i(x)L - 2H lihlfijl (x)L 
+ H 0Xltl (x)L 2 + H m {x)L 2 + H ltlfl>l {x)L 2 - H w (x)L 2 - ^n 4 H 1 (x)L 

+ l7r 2 // ,i,i(x)L + i7T 2 //i, ,i(x)L + 7T 2 H 1:1A (x)L + y^L 2 - ^n 2 H 0A (x)L 2 

- ^n 2 H hl (x)L 2 + ^tt 2 //i(x)L 3 - ^ ,i,i(x)i 3 - ^i, ,i(x)L 3 + |ffi,i,i(x)L s 

- ^ 4 + -^o,i(*)^ - ^i,i(^)L 4 + ^i(x)L 5 - 1l 6 + ^ 4 #i,i(*) 
1111 

+ g7T 2 #0,0,0,l(^) + -7T 2 i7 ,0,l,l(^) + g?T 2 #0,1,0,1 (x) + -7T 2 ifi, ,0,l (x) ~ 7r 2 #i,i,i,i (x) 

5 14 26 14 

- -vr 2 iJi( a :)C3 - yiWaOCa + 2ffi,i,i(x)C 3 + y^Cs - y #o,i(^C 3 

4 4 1 40 72 1 

- 2H hl (x)L( 3 - -#iOr)L 2 C 3 + g^Cs + -yCf + 2^i(x)C 5 - y^5 + ^tt 6 . (B3) 



2. Two- loop integrals 



The two-loop planar scalar double-box integral depicted in fig. 7 is 
li 2 \s,t) = (-ie"*-** 2 ) 2 J 



d d pd d q 



p 2 (p — ki) 2 (p — ki — k 2 ) 2 (p + g) 2 g 2 (q — k^) 2 (q — k 3 — k^) 2 

(B4) 

This integral was first evaluated in ref. [14] through 0(e°), as required in NNLO calcula- 
tions. Here we need the integral through 0(e 2 ). The calculation performed in ref. [14] was 
not optimal because the starting point was a fivefold MB representation. On the other 
hand, it is possible to derive an appropriate fourfold representation, as was demonstrated in 
ref. [69] (see also Chapter 4 of ref. [48]). The corresponding evaluation can be generalized 
straightforwardly to obtain the next two orders of the expansion in e. Let us stress that this 
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evaluation is much simpler than the evaluation of the triple boxes up to e c 
Our result through C(e 2 ) is, 



1 Cj (x, L) 



where x = —t/s, L = ln(s/£), and 

5 

c 4 = -4 , c 3 = -5L , c 2 = -2L 2 + -vr 2 , 

2 11 6 5 

Cl = 4 [-L# ,i(a:) " #o,o,i(a:)] — 2 (L 2 + vr 2 ) if^x) + -L 3 + —Lit 2 + -( 3 , 
c = 4 [lliZ ,o ) o,i( a O + # ,o,i,i(a:) + #0,1,0,1 ( x ) + #1,0,0,1(2:)] 
+4L [6^0,0,1 (a:) + ^i.iW + #1,0,1(2:)] + 2 ^ 2 [#o,i(a:) + ffyW] 
+ 2 -7i 2 [10H 0A (x) + 3#i,i(:r)] + ^H^x) [-4L 3 - 5Ltt 2 - 6(3] 

c_i = -4 [28^0,0,0,0,1 + 29^0,0,0,1,1 ( a: ) + 24# ,o,i,o,i(a;) + #0,0,1,1,1 (x) 
+19# ,i,o,o,i(a;) + #0,1,0,1,1(3;) + #0,1,1,0,1(3:) + 14-ffi,o,o,o,i (x) 
+#1,0,0,1,1 (x) + #1,0,1,0,1(3:) + #1,1,0,0,1(3:)] 

-4L [18^0,0,1,1(3:) + 13^0,1,0,1(3:) + #0,1,1,1(3:) + 8H 1>0fitl (x) + #1,0,1,1(2;) + #i,i, 
+2L 2 [12# ,o,i(x) - 7^0,1,1 (x) - 2#i, ,i(:r) - #1,1,1(2:)] 
^vr 2 [#0,0,1(2:) - 28# ,i,i(x) - 13#i, ,i(2:) - 3#i,i,i(x)] 

j|L 3 [2# ,i(x) + #i,i(x)] + ^Ltt 2 [18# ,i(x) + 5#i,i(x)] 

+72C 3 #o,i(2:) - 4 [36L 4 + 78LV + 17tt 4 ] #i(x) - 4( 3 [-7LH x {x) - H hl (x)\ 

to 

14 r5 13 2r3 46 A r2 211 4r 73 2a 383 a 

H L H 7r 2 L 3 H (3^ H tt 4 L vr 2 C 3 H Cs , 

1 R : 2 3 s3 1 90 fi S3 Fi S5 ' 



+ 
+ 



15 3 3 J 120 6 J 5 J ' 

c„ 2 = 4 [68# ,o,o,o,o,i(2;) + 76# , ,o,o,i,i(3;) + 66^/0,0,0,1,0,1 (») + 65# , ,o,i,i,i(3;) 

+56# ,o,i,o,o,i(2;) + eOf/'o.o.i.o.i.i (») + 48# , ,i,i,o,i(a:) + #0,0,1,1,1,1 (a:) 

+46# ,i,o,o,o,i(a:) + 55# ,i, ,o,i,i(a:) + 43# ,i, ,i,o,i(a:) + #0,1,0,1,1,1 (a:) 

+31# ,i,i,o,o,i(a:) + #0,1,1,0,1,1 (a:) + #0,1,1,1,0,1(2:) + 36#i, ,o,o,o,i(a:) 

+50#i )0 , ,o,i,i(a:) + 38#i, , ,i,o,i(a:) + #1,0,0,1,1,1 (a:) + 26#i, ,i, ,o,i(a:) 

+#1 ,0,1,0,1,1(2:) + #1,0,1,1,0,1 (a:) + 14#i, 1,0,0,0,1 (x) + #1 ,1,0,0,1, 1(3;) 

+#1,1,0,1,0,1(2:) + #1,1,1,0,0,1(2:)] + 4L [42# ,o,i,i,i(x) + 37# ,i,o,i,i(a:) 
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+25# ,i,i,o,i(z) + #o,i,i,i,i(z) + 32# 1 , , ,i,i(x) + 20^,0,1,0,1 (x) + H lfitlxl (x) 
+8#i,i, ,o,i(z) + #i,i,o,i,i(z) + ^1,1,1,0,1^)] " 2L2 [36# ,o,i,i(z) + 26# ,i,o,i(z) 

-19^0,1,1,1^) + MH lt0fitl (x) - 14^,0,1,1 (x) - 2H hlA1 (x) - #i,i,i,i(x)] 
2 

--7T 2 [17^0,0,0,1^) + H ,o,i,i(x) + 6# ,i,o,i(x) - 64iZ o ,i,i,i( a 
+11^1,0,0,1^) - 49ifi i0) i 1 i(aO - 13//i,i, ,i(x) - 3#i,i,i,i(x)] 
+^L 3 [6tf ,o,i(z) - 8fl- 0l i,i(x) - 3^ li0) i(x) - H lxl (x)] 

~n 2 L [6H ,o,i(x) + 54iZ O) i,i( a + 29# 1 , ,i(z) + Stf^,^)] 

-^Cs [40tf ,o,i(z) + 90fZ- 0) i,i(x) + 75H lfiA (x) + Stf^,^)] + h 4 [7H 0<1 (x) + 3/^Or)] 
+^ 2 L 2 [33H 0tl (x) + 13^,!^)] + ^vr 4 [129H 0>1 (x) + SS^,^)] 
+^[59^ 0l i(a:)-21ff 1)1 ( a ;)] 

+ — [-48L 5 - 160LV - 55Ltt 4 + 1140L 2 C 3 + 240tt 2 C 3 - 720C 5 ] H x (x) 
45 

+ + 1^ + ^ + _ 

Through O(e ) this corresponds to the results of ref. [14]. 

It is also possible to derive differential equations obeyed by the planar two-loop box 
integral [21]. The differential equations couple I^\s,t) to a second master two-loop box 
integral. In ref. [21] these results were used to obtain the second integral, and to check 
I^\s,t) through order e°. We have used the same differential equations to check the re- 
sult (B6) through the required order, e 2 , up to a constant. The order e 2 constant was checked 
numerically. 



3. Three-loop integrals 

The three-loop ladder integral depicted in fig. 4(a) and defined in eq. (3.1) has been 
evaluated in ref. [23], with the result, 

v ; 3=0 
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where x = —t/s, L = ln(s/t), and 

16 5 r 3 2 

ce = y, c 5 = --L, c 4 = --7T , 

3 11 131 

c 3 = 3(^0,0,1^) + ^o,i(^)) + 2 ( L2 + O^iO*) - i2 rt - — ^ 3 > 

c 2 = -3 (17if ,o,o,i(^) + #0,0,1,1(3:) + H 0>lt o tl (x) + #1,0,0,1(3:)) 
—L (37H 0fi>1 (x) + 3ff ,i,iW + 3^i )0l i(x)) - ^(L 2 + n 2 )H 1A (x) 

- (f L 2 + 8tt 2 ) // ,i(x) - Ql 3 + tt 2 L - 3C 3 ) H x {x) + f(,L - ^ , 
ci = 3 (81^0,0,0,0,1^) + 41f/ ,o,o,i,i( a; ) + 37i7o,o,i,o,i(z) + #0,0,1,1,1 (a) 
+33ifo,i,o,o,i(a0 + #0,1,0,1,1(2:) + #0,1,1,0,1(2;) + 29^1,0,0,0,1 C^) 
+#1 ,0,0,1, + #1 ,0,1,0,1(3:) + #1,1,0,0,1(2;)) + L (177ifo,o,o,i(^) + 85f/o,o,i,i(^) 
+73if 0) i,o ) i(aO + 3^0,1,1,1^) + ^H 1>0fi>1 (x) + 3ifi )Ol i,i(a0 + 3H 1>lfi>1 (x)) 
+ + ^f- 2 ) H 0fitl (x) + (yL 2 + 20vt 2 ) // ,i,i(x) 

+ (|L 2 + 14vr 2 ) H^x) + \ (L 2 + vr 2 ) ^(x) 

+ (f # 3 + §tt 2 L - 96C 3 ) #0,1 (*) + (|l 3 + vr 2 L - 3C 3 ) #1,1 (s) 

/9 r4 25 2r2 rf1> r 13 A TT , N 503 4r 73 2> 301, 

c = - (951# ,o,o,o,o,iOe) + 819i? o ,o 1 o,o,i,i( a + 699-H Ol o,o 1 i,o,i( a + 195if ,o,o,i,i,i(^) 
+547^0,0,1,0,0,1(3?) + 231iJo,o,i,o,i,i(aO + 159^0,0,1,1,0,1^) + 3-H 0l o ) i,i,i,i( a: ) 
+363^0,1,0,0,0,1 (x) + 267if ,i,o,o,i,i(^) + 195if ,i,o,i,o,i(^) + 3if ,i,o,i,i,i(^) 
+123// ,i, 1,0,0,1 + 3 #o,i, 1,0,1, i(x) + 3-fZo,i,i,i,o,i(x) + 147ifi, , ,o,o,i(x) 
+303^1,0,0,0,1,1 + 231ifi,o,o,i,o,i(aO + 3if 1,0,0,1,1,1 (a?) + 159i7i, ,i,o,o,i(^) 
+3^1,0,1,0,1,1(3:) + 3ffi )0l i,i,o,i(a:) + 87ifi, 1,0,0,0,1 (x) + 3i7i, 1,0,0,1, i(x) 
+3i^i,i,o,i,o,i(3:) + 3ifi, 1,1,0,0,1(3:)) 

-L (729H omi (x) + 537^0,0,0,1,1(2:) + 445i/o,o,i,o,i(3:) + 133^0,0,1,1,1^) 

+321/7 ,i,o,o,i(3:) + 169/7 ,i,o,i,i(3:) + 97/7 ,i,i,o,i(3:) + 3iJo,i,i,i,i(3:) 

+165//i, ,o,o,i(3:) + 205//i, ,o,i,i(3:) + 133//i, ,i, ,i(3:) + 3//i, 0,1,1,1(3:) 

+6I//1, 1,0,0,1(3:) + 3//i, 1,0,1, 1(3;) + 3Z/i, 1,1,0,1(3:)) 
/531 r2 89 2 \ r , , /311 r2 619 2 \ 
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- (^ 2 + ^- 2 ) %mW - (y^ 2 + 32vr 2 ) ff„, ul W 

- - ^- 2 ) #i,o,o,i(*) - + 50^) #i,o,i,i(*) 



- (f L 2 + 14vr 2 ) ff 1>1>0>1 (x) - \ (L 2 + vr 2 ) if^s) 

- (^ 3 + ^- 2 L - 455C 3 ) - (y^ 3 + ^tt 2 L - 120C 3 ) #0,1,1' 

- (f ^ 3 + - 156C 3 ) #i,o,i(*) - Ql 3 + - 8(3) #i,i,i(*) 

/69 r4 101 2r2 r 559 ,\ , v 

(^y^ 4 + — ^ 2 ^ 2 - 291C 3 £ + ffo,i(aO 

-gL 4 + | 7 r 2 L 2 -58C 3 L + ^i/ 1 , 1 ( a ;) 

27 r5 25 2r3 183 2 131 4 r 37 2a r ^ \ TT . . 
40 L + ~8~ - — CaL + ~60~ " 12^ + 5?Cs ) H ^ X) 

223 2^ n r 167 ^2 624607 6 

ir 2 C3+149 (5 L+^Cf- —n>. 



The result for the second triple box, defined in eq. (3.2) and shown in fig. 4(b), 

4 (-S) 1+3e t 2 ^ 6* ' 

where 



16 13 T 1 r2 19 2 

Cfi = — , C5 = — L , C4 = — L 71 , 

6 9 ' 6 ' 2 12 ' 

c 3 = I [#0,0,1 (*) + £#0,1 (a:)] + | [£ 2 + ^r 2 ] #i(x) 

-1l3-^ 2 -^ C 3, 
12 72 18 

c 2 = - [ll#o,o,o,i(aO - 5 #0,0,1,1(20 ~ 5^0,1,0,1 (x) - 5H lfifi ,i(x)] 

+U [UH 0fiA (x) - 5# ,i,i(z) - 5ifi i0) i(a;)] + \l 2 [17H 0i1 (x) - 5H hl (x)} 

+^H , 1 (x) - ^H 1A (x) + ^L 3 Hi(x) + ^Ln'H^x) 

41 r, 5 , . l r4 l r2 2 2429 4 
-T LC * + 2 H ^-3 L -4 LV + 6480^' 

ci = - [-55^0,0,0,0,1 ~~ 59if 0,0,0,1,1 (z) - 31^0,0,1,0,1 W + 5# o ,o,i,i,i(a0 
-3# ,1,0,0,1 w + 5# ,1,0,1,1 (a;) + 5# ,1,1,0,1(2;) + 25/^1,0,0,0,1 (#) 
+5-ffi,o,o,i,i + 5 #i,o,i,o,i (a;) + 5#i,i, ,o,i(z)] 
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+h [22H 0tOfitl (x) - A6H 0fiX1 (x) - 18H 0ilfiil (x) + 5H 0jl;ljl (x) 
+10^1,0,0,1 (x) + 5//i, ,i,i(x) + 5//i,i, ,i(x)] 
+^L 2 [64// ,o,i(^) - 33# ,i,i(z) - 5//i, ,i(x) + 5ifi,i,i(a;)] 

^tt 2 [25^0,0,1 (^) " 128^ 0l i,i(^) + 40#i, ,i(z) + 30H 1AA (x)] 
±L 3 [71H 0A (x) - 20H 1A (x)} 

■^Lk 2 [153^0,1^) " 50^i,i(x)] + l - [8H 0A (x) - 5H ltl (x)] ( 3 
43 71 ^ ^ 7 

48 v ; 48 v ; 144 v ; 2 v ;s 48 

227 . 2 13 r2> 10913 r 4 3257 2a 889 , 

H L V H L 2 C 3 H Lit 4 H 7T 2 C 3 C5 , 

^144 4 S3 8640 216 S3 10 s5 ' 



+ 

1 

+ 
+ 



144 4 Sd 8640 216 s 10 s ' 

= X - [379^0,0,0,0,0,1 (^) + 343// ,o,o,o,i,i(^) + 419// ,o,o,i,o,i(^) + 347// ,o,o,i,i,i(^) 
+355i? ,o,i,o,o,i(^) + 175i? ,o,i,o,i,i(^) + 223i? ,o,i,i,o,i(^) - 5 #0,0,1,1,1,1(2;) 
+151i? ,i,o,o,o,i(^) + 3^0,1,0,0,1,1 (a?) + 51i? ,i,o,i,o,i(^) - 5i? ,i,o,i,i,i(^) 
+99i? ,i,i,o,o,i(a:) - 5 #0,1,1,0,1,1(2;) - 5# ,i, 1,1,0,1 ( x ) ~ 193#i,o,o,o,o,i(a;) 
-169iJi, ,o,o,i,i(z) - 121-H"i, ,o,i,o,i(^) - 5-f/"i, ,o,i,i,i(^) - 73ifi, ,i, ,o,i(^) - 5 #1,0,1,0,1,1(2:) 
~~ 5-£/i,o,i, 1,0,1 (^) — 25i?i, 1,0,0,0,1 (x) ~~ 5i?i, 1,0,0,1,1 (x) — 5#i, 1,0,1,0,1(2;) — 5#i, 1, 1,0,0,1 (2:)] 
+h [98H omi (x) - 22^0,0,0,1,1(2:) + 98# ,o,i,o,i(z) + 238// ,o,i,i,i(^) + 78iJ ,i,o,o,i(^) 
+66// ,i,o,i,i(^) + 114//o,i,i,o,iW - 5# ,i,i,i,i(z) - 82^1,0,0,0,1(3;) - 106^1,0,0,1,1(0:) 
-58i?i, ,i,o,i(a:) - 5ifi j0 ,i,i,i(a;) - 10#i,i, ,o,i(2;) - 5if 1,1,0,1,1(2:) - 5#i,i,i,o,i(a;)] 
+h 2 [12AH 0fiAl (x) - 208^o,o,i,i (x) - 44# , 1 , , 1 (x) + 129/J ,i,i,i(^) 
-20^1,0,0,1 (x) - 43ifi )0l i,i(a;) + 5#i, 1 , ,i(z) - 5H hW {x)] 
+^tt 2 [183iZ 0) o,o ) i(aO - 121/J ,o,i,i(^) + 375# ,i,o,i(2;) + 704/J ,i,i,i(^) 
+31^1,0,0,1(0;) - 328//i,o,i,i(x) - 40 J ffi,i )0l i(x) - 30ifi,i,i,i(a;)] 
+^L 3 [260// ,o,i(o:) - 215^o,i,i(a:) - 7H lfitl (x) + 20if 1 , 1)1 (a;)] 
+ ^Ltt 2 [326H , ,i(x) - 633if o ,i,i(a0 + 127#i, ,i(2;) + 50^i,i,i(x)] 
l - [-3LH 0A {x) - 5LH 1A (x) + 165H 0fiA {x) + 104if 0l i,i(aO - 68#i, ,i(:r) - 5ifi,i,i(a;)] Ca 
^L 4 [309fZ- Ol i(x) - 43ff 1)1 (x)] + ^LV [725if 0) i(aO - 71#i, 1 (o:)] 
A [1848H 0A (x)+25H 1A (x)] 



+ 
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+ ^ Hl ( X ) + H L »^ l(l ) + ^Lt^s) + f ^3 + ^^3 

-2L 2 ^(-)C3 - ^Wi(*)Ca + ^iMCs + ^ + ||lV 

3163 24 1054 1391417 6 197 2 " 

+ 2880 L ff " ^ LCs + 3265920" + ~6~^ 3 " (B10) 
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